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We study single-ion and exchange anisotropy effects in equal-spin s\ tetramer single molecule 
magnets exhibiting Td, D^h, Did, Cah, Ca v , or 54 ionic point group symmetry. We first write the 
group-invariant quadratic single-ion and symmetric anisotropic exchange Hamiltonians in the ap- 
propriate local coordinates. We then rewrite these local Hamiltonians in the molecular or laboratory 
representation, along with the group-invariant Dzyaloshinskii-Moriya (DM), and isotropic Heisen- 
berg, biquadratic, and three-center quartic Hamiltonians. Using our exact, compact forms for the 
single-ion spin matrix elements, we evaluate the eigenstate energies analytically to first order in the 
microscopic anisotropy interactions, corresponding to the strong exchange limit, and provide tables 
of simple formulas for the energies of the lowest four eigenstate manifolds of ferromagnetic (FM) 
and anitiferromagnetic (AFM) tetramers with arbitrary si. For AFM tetramers, we illustrate the 
first-order level-crossing inductions for si = 1/2,1,3/2, and obtain a preliminary estimate of the 
microscopic parameters in a N14 from a fit to magnetization data. Accurate analytic expressions for 
the thermodynamics, electron paramagnetic resonance absorption and inelastic neutron scattering 
cross-section are given, allowing for a determination of three of the microscopic anisotropy interac- 
tions from the second excited state manifold of FM tetramers. We also predict that tetramers with 
symmetries S4 and D^d should exhibit both DM interactions and multiferroic states, and illustrate 
our predictions for s\ = 1/2, 1. 

PACS numbers: 75.75.+a, 75.50.Xx, 73.22.Lp, 75.30.Gw, 75.10.Jm 



I. INTRODUCTION 

Single molecule magnets (SMM's) have been a topic 
of great interest for more than a decade, [lj because of 
their potential uses in quantum computing and/or mag- 
netic storage, Q which are possible due to magnetic quan- 
tum tunneling (MQT) and entangled states. In fits to a 
wealth of data, the Hamiltonian within an SMM cluster 
was assumed to be the Heisenberg exchange interaction 
plus weaker total (global, or giant) spin anisotropy inter- 
actions, with a fixed overall total spin quantum number 
s.[H MQT and entanglement were only studied in this 
simple model. 

The simplest SMM clusters are dimers.@, 0, Q 
Surprisingly, two antiferromagnetic dimers, 
an Fe2, [Fe(salen)Cl]2, where salen is N,N'- 
ethylenebis(salicylideneiminato), and a M2, 
Na2Ni2 (0204)3(1120)2, appear to have substantial 
single-ion anisotropy without any appreciable total 
spin anisotropy. @, g, Q The presence of single-ion or 
exchange anisotropy actually precludes the total spin 
s from being a good quantum number. 0, [B| Although 
the most common SMM clusters have ferromagnetic 
(FM) intramolecular interactions and contain n > 8 
magnetic ions,[H, Q a number of intermediate-sized FM 
SMM clusters with n = 4 and rather simple molecular 
structures were recently studied. Fits to electron 
paramagnetic resonance (EPR) N14 data assuming a 
fixed s were also problematic, suggesting single-ion or 
exchange anisotropy in that tetramer, as well, jiol 11 1 



splittings. [TJ, 
where hmp 



The Cu 4 tetramer Cu 4 OCl 6 (TPPO) 4 , where TPPO is 
triphcnylphosphinc oxide, has four spin 1/2 ions on 
the corners of a regular tetrahedron, with an s = 2 
ground state and approximate Td symmetry. 12|, Il3l . 1 14 1 
In this case, there are no single-ion anisotropy effects, 
but anisotropic symmetric exchange interactions were 
thought to be responsible for the zero-field energy 
The Co 4 , Co 4 (hmp)4(MeOH) 4 Cl4, 
is hydroxymethylpyridyl, and Cr4, 
[Cr4S(02CCH3) 8 (H20)4](N0 3 )2-H 2 0, compounds 
have s = 6 ground states with spin 3/2 ions on the cor- 
ners of tetrahedrons . (Til. [l7| Those compounds have S4 
and approximate Did symmetry, respectively. 1(1 17 1 A 
number of high symmetry s — 4 ground state Nu struc- 
tures with spin 1 ions were reported. [H, [l9|, [2Q 1 , [U |22[ 
Two of these, [Ni(hmp)(ROH)Cl] 4 , where R is an alkyl 
group, such as methyl, ethyl, or 3,3-dimethyl-l-butyl 

form tetramers 
Two others, 

Ni 4 (ROH)Ij4, where R is methyl or ethyl and H 2 L 
is salicylidene-2-ethanolamine, had approximate S4 
symmetry, although the precise symmetry was only 
Ci.[3 Several planar Mn4 compounds with the Mn +3 
spin 2 ions on the corners of squares were made, with 
overall s = 8 tetramer ground states. [23] Although two 
of these complexes had only approximate S4 symmetry, 
one of these complexes, Mn4Ci4(L')4, where H2L' 
is 4-i-butyl-salicylidene-2-ethanolamine, had perfect 
S4 symmetry. [23j Inelastic neutron scattering (INS) 
experiments provided strong evidence for single-ion 



and hmp is 2-hydroxymethylpyridyl 
with precise S4 group symmetry. pll |2 



2 



anisotropy in a C04 
symmetry. 



17, 1 



and a M4 with approximate S4 



We note that ab initio calculations of the intramolec- 
ular spin-spin interactions in SMM clusters have not yet 
been always successful in calculating even the strongest, 
intramolecular isotropic Heisenberg interactions accu- 
rately, and have been incapable of calculating any of the 
local anisotropi c sp in-spin interactions within an SMM 
cluster. [U [23, 0| Even to obtain the Heisenberg in- 
teractions accurately, it seems one needs to extend the 
local spin-density approximation (LSDA) to include on- 
site repulsions with strength U (the LSDA+U model), 
which would have to be introduced phcnomcnologically 
to fit the lowest two energy level manifolds in zero applied 
magnetic field. 0, [27|, S [H, [1| We therefore define a 
microscopic model to be a model constructed in terms of 
the individual spins and from the local interactions be- 
tween them, with parameters describing the strengths of 
the various types of local spin-spin interactions and in- 
teractions between the local spins and the magnetic field. 
This is distinct from a model constructed solely from the 
anisotropies of the total spin of an SMM cluster, which we 
denote as a phenomenological model. Our definition of a 
microscopic model is analogous to the standard model of 
the interactions of quarks and gluons within a hadron. 

Recently there have been microscopic treatments of 
dimers,[3,li| trimers, and tetramers, including Zeeman ex- 
tensor anisotropy, single-ion anisotropy, and anisotropic 
exchange interactions. [3 lj Most of those treatments and 
their recent extensions to more general systems expressed 
the single-spin matrix elements only in terms of Wigner 



3j, 6j, and 9j symbols. [3 ll |32| While such treatments 



are very helpful in fitting experimental data, more com- 
pact analytic forms are desirable to study microscopic 
models of FM SMM clusters in which the MQT and 
entanglement issues crucial for quantum computing can 
be understood. We constructed the quadratic single-ion 
and anisotropic near-neighbor (NN) and next-nearest- 
neighbor (NNN) exchange SMM cluster Hamiltonians 
from the respective local axial and azimuthal vector 
groups for equal-spin tetramer SMM clusters with point 
group symmetries g = T d , D ih , D 2 d, C ih , C iv , and £4, 
and found compact analytic expressions for the single- 
spin matrix elements of four general spins. Each local 
vector group generates site-dependent molecular single- 
ion and exchange anisotropy. We then show that for 
and S4 symmetries, the antisymmetric exchange inter- 
actions lead to non-vanishing spin currents that may be 
accompanied by electric polarizations, leading to multi- 
ferroic effects. We evaluate the magnetization, specific 
heat, EPR and INS transitions in the Hartree approxi- 
mation, and provide a procedure for extracting three of 
the effective site-independent microscopic parameters us- 
ing EPR. We also show analytically how to include the 
effects of weak biquadratic exchange. 

An outline of the paper is as follows. In Sec. II, we 
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FIG. 1: T d (left) and D 4h 

origin. Arrows: local axial z„ d (left), azimuthal x 
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ionic plane. 
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FIG. 2: D 2d (left) and S4 (right) ion sites (filled). Circle: ori- 
gin. Arrows: local axial single-ion vectors. The g — D2d,Si 
axial vectors zf make the angles 8f with the z axis, and the 
S4 axial vector zf 4 also makes the angle a\ with the x axis, 
where cos at = sin 6f 4 cos cjf>f 4 . 



discuss the six structures and the general quadratic spin 
Hamiltonian. In Sec. Ill, we write the single-ion and 
symmetric anisotropic exchange Hamiltonians in terms 
of the local coordinates, and the antisymmetric exchange 
Hamiltonian in the molecular coordinates. In Sec. IV, 
we impose the operations of the six group symmetries, 
and discuss the effects of antisymmetric anisotropic ex- 
change interactions and the related electric polarizations 
in lower symmetry systems. In Sec. V, the resulting 
group-symmetric Hamiltonians are written in the molec- 
ular representation, and the isotropic biquadratic ex- 
change interactions are introduced. Section VI contains 
the eigenstates of the full Hamiltonian to first order in the 
anisotropy and NN biquadratic exchange interactions. 
These eiqenstates are used to obtain the level-crossing 
inductions for AFM tetramers, and particular examples 
with s\ = 1/2,1,3/2 are presented. In Sec. VI, we 
also evaluate quantitatively some effects of antisymmet- 
ric anisotropic exchange and provide our related predic- 
tions for multiferroic behavior. In Sec. VII, the self- 
consistent Hartree approximation (or strong-exchange 
limit) is used to provide simple but accurate results for 
the thermodynamics, EPR resonant inductions, and INS 
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g = Did or 5*4, c/a > 1, approximately as in a C04, 17 1 
or c/s < 1, as in some Mn4 and a M4. TjJ 23| For com- 
parison with the planar symmetries g = C^D^, and 
D± v , we assume for g — D 2 d,S4 that c/a < 1, so that 
there are four NN sites and two NNN sites. For each g, 
x, y, z are the molecular (or laboratory) unit coordinate 
axis vectors. 

The most general Hamiltonian quadratic in the four 
spin operators S n may be written for group g as 



FIG. 3: C ih (left) and C iv (right) ion sites (filled). Circle: 
origin. Arrows: local azimuthal Xn ih (left) axial z^ 4 " (right) 
single- ion vectors. The axial vectors z n ih = z. The d v axial 
vectors each make the angle 61 = n/2 — 6[ with the z axis. 
The dotted arrows (equivalent to the D^h azimuthal single-ion 
vectors x n ih ) are their projections in the xy plane. 



cross-sections, and describe how EPR experiments in the 
excited states of FM tetramers can provide a measure of 
some of the microscopic anisotropy interations strengths. 
Finally, in Sec. VIII, we discuss the significance of our re- 
sults, and provide a preliminary fit to magnetization data 
on an AFM N14 tetramer, and in Sec. IX, we present our 
conclusions. 



II. STRUCTURES AND BARE HAMILTONIAN 

For SMM clusters with ionic site point groups g = 
Td,D4h, we assume the four equal-spin si ions sit on 
opposite corners of a cube or square of side a centered 
at the origin, as pictured in Fig. 1. For clusters with 
g = D 2 d, S4, we take the ions to sit on opposite corners 
of a tetragonal prism with sides (a, a, c) centered at the 
origin, as in Fig. 2. The ions for g = C^h, Civ also sit on 
the corners of a square of side a centered at the origin, 
as pictured in Fig. 3, but the ligand groups have differ- 
ent symmetries than for the simpler D±h case pictured 
in Fig. l.[33| In each case, we take the origin to be at 
the geometric center, so that J2 n =i r ™ = 0, wnere the 
relative ion site vectors are 



r n = 



a 

V2 



. f{2n— 1)tt\ „ / 
smy Jx + cosy 



(2n- 1)tt 



-(-!)"£. 



(1) 



Tetrahedrons with g — Tj, c/a — 1, approximately as 
in Cu4,[T2| are a four-spin example of the equivalent- 
neighbor model. [34:] In squares with g = D^h, C^h, or 
Ci V , c = 0. The high symmetry is approximately ex- 
hibited by the square Nd4 compound, Nd4(OR)i2, where 
R is 2,2-dimethyl-l-propyl, in which the Nd +3 ions have 
equal total angular momentum j = 9/2.[35|, H(| We note 
that the Mn4 clusters with approximate or exact S4 sym- 
metry also have c = 0.[23] In tetragonal prisms with 



«' = -^E B - 9 9 n -S n + S " ■ D n,n> ' ^V, (2) 



where /is is the Bohr magneton and B = 
B(sm9cos(f>,sm9sm(j),cos9) is the magnetic induc- 
tion at an arbitrary direction (0, (f>) relative to the 
molecular (or cluster) coordinates (x. v. z).(3ll. I37I] 

For simplicity, we take *g g n to be diagonal, isotropic, 
and site-independent, so that the Zeeman interaction 
may be written in terms of a single gyromagnetic ra- 
tio 7 ~ 2fis- Thus in the following, g only refers to 
the molecular group. We separate D 9 n , into its sym- 
metric and antisymmetric parts, D 9 n n , = D 9 ^ s n , + D 9 ^ a n ,, 
respectively. For n' — n, the single-ion D 9 n n is neces- 
sarily symmetric, so D-fc^ = 0. For each g, the four 
Dn S n contain the local single-ion structural information, 

and the six distinct symmetric D 9 ' s , contain the local 
symmetric exchange structural information, which lead 
to the isotropic, or Heisenberg, exchange interactions, 
and the remaining symmetric anisotropic exchange in- 
teractions. The six distinct antisymmetric D 9 /^, con- 
tain additional local structural information which lead 
to the Dzyaloshinskii-Moriya (DM) interactions. [38l [39j] 
Physically, the symmetric anisotropic exchange interac- 
tions also contain the intramolecular dipolc-dipolc inter- 
actions, which can be even larger in magnitude than the 
terms originating from actual anisotropic exchange. 37. 

E3 

As is well known, each of the symmetric rank-three 
tensors (or matrices) D 9 n n , can be diagonalized by three 
rotations: a rotation by the angle (j) 9 n n , about the molec- 
ular z axis, then a rotation by the angle Q 9 n , about the 
rotated x axis, followed by a rotation by the angle i\> 9 n 

about the rotated z axis, 

^3 



41] This necessarily leads to 
the three principal axes x n n , , y nn i, and z n n , . For the 
single-ion axes with ml = n, we denote these principal 
axes to be x n , y„, and z n , respectively, which are writ- 
ten explicitly in Sec. III. The non-vanishing matrix ele- 
ments in these locally-diagonalized symmetric matrix co- 
ordinates are D 9 ji, „„, D 9 n s n , yy and D^ n , >zz . Since the 



structural information in each of the D 9 / S n , depends upon 
the local environment, in the absence of molecular group 
g symmetry, each of these angles would in principle be 



4 



different from one another. 

Although an antisymmetric exchange matrix D^Zi can 
generally be diagonalized by a unitary transformation, 
it contains at most three independent, real parameters, 
which can be incorporated into the components of a 
three- vector, d a nn ,, with an effective spin-spin interac- 
tion of the form d 9 n n , ■ (S„ x S„> ) , [H, El which is easiest 
to write in the molecular representation. 

For the six high-symmetry groups under study, we an- 
alyze the effects of molecular group symmetry upon the 
single-ion and anisotropic exchange parts of 7i 9 . The 
group symmetries further restrict the number of inde- 
pendent parameters. 

In the absence of any anisotropy interactions, the bare 
Hamiltonian TLq is given by the Zeeman and Heisenberg 
interactions, 



basis elements are the vectors that diagonalize the single 

ion matrix from D 9 n n to D 9 L „.[4l[ Since we employ these 
vectors repeatedly, we write them here for simplicity of 
presentation. The diagonalized vector set elements may 
be written in the molecular (x,y,z) representation as 



■En 



%9 



Z,„ = 



cos #2 sin § 9 n sin tp^ 



cos<^cosV>k •„ 
sin <p 9 n cos ip 9 + cos 9^ cos (p 9 n sin ip> 
sin 9^ sin ips 



— cos <p 9 n sin ip9 — cos 9^ sin <p 9 n 

sin cp 9 x sin tp^ + cos 9^ cos § 9 n cos ip^ 



sin fib sin tp 9 ^ + cos 9^ cc 
sin 9 9 cos ipn 



< 9 COS 1pn 



sin 6*9 sin 09 
- sin 9 9 cos (p 9 n 
cos 6 9 n 



(7) 




H% = -jB-S- J' g (S 1 -S 3 + S 2 -S 4 ) 

-J g (Si ■ S 2 + S 2 - S 3 + S 3 - S 4 + S 4 - S x ), (3) 



which can be rewritten as 
J, 



n 



{J' q - Jg) . 



fS'-^B-S- ^-2-^(S?3 + (4) 



where Si 3 = Si+S 3 , S->4 = S2+S4, and S = S13+S24 is 
the total spin operator, [36( and we dropped an irrelevant, 
overall constant. In Eq. Q, 



J r d - Jr d , 

Jg 7^ Jg 



(5) 
(6) 



for g = D2d, S4, D4h, C4h, and Cav. In terms of the diag- 
onalized matrix elements, —2J 



■^1,2, xx 



D 1,2, VV and 



— 2J' g = Dl'l + Df'l yy , for instance. For our c/a < 1 
convention, — J g and — J' g are the NN and NNN Heisen- 
berg interactions for g = D 2 d, S4, C4 Vl C^h and D4h- 



III. THE SINGLE-ION AND ANISOTROPIC 
EXCHANGE HAMILTONIANS 

To take account of the molecular group g symme- 
tries, it is useful to write the single-ion and symmetric 
anisotropic exchange interactions in terms of the local 
coordinates. In this section, we write the local Hamil- 
tonian for these interactions, and the molecular Hamil- 
tonian for the antisymmetric exchange interactions. In 
Sec. IV, we then impose the group symmetries on these 
interactions for C^, -D^, C4„, S4, Did-, and Td molecular 
group symmetries, respectively. 



A. Local single-ion Hamiltonian 



For the single-ion anisotropy, we define the local vector 
basis for the nth site to be {x n , y n , z n } for each g. These 



which satisfy x 3 n x y n = z n . We then write the most 
general quadratic single-ion anisotropy interaction as 

4 



+J?, n [(S n -x. 9 n f-(S n .ti 9 n ) 2 }), (10) 

in terms of the site-dependent axial and azimuthal in- 
teractions JanjJiju analogous in notation to that for 
homoionic dimers.0, 0] In terms of the diagonalized ma- 
trix elements, -J£„ = D^ n>zz - (D 9 n ^ xx + D 9 n ^ yy )/2 
and -Jg = 



n,n,xx 



D 9 ' s )/2 

n,n,yy ) I 



Local symmetric anisotropic exchange 
Hamiltonian 



In addition to the single-ion interactions, the other 
microscopic anisotropic interactions are the anisotropic 
exchange interactions, which include the intracluster 
dipole-dipole interactions. The intercluster dipole- 
dipole interactions can lead to low-T hysteresis in the 
phenomenological total spin model. (4^| but in the micro- 
scopic individual spin model, are generally much weaker 
than the intracluster ones due to the larger distances in- 
volved. Hence, we neglect those and all other intercluster 
interactions, such as those mediated by phonons. As for 
the single-ion interactions, we first construct the sym- 
metric anisotropic exchange Hamiltonian TL 9 ae in the lo- 
cal group coordinates. In this case, there are distinct 
local vector sets for the NN and NNN exchange inter- 
actions. Diagonalization of the symmetric anisotropic 

exchange matrix D„ a n , leads to D 9 ^ s n , and the vector ba- 
sis {Xn,n'>li 9 n,n'i2 9 n,n'}i S iven by Eqs. ©-© with the 
subscript n replaced by n,n'. 

The local symmetric anisotropic exchange Hamiltonian 



5 



W^'g is then generally given by 

2 6-2g 



q— 1 n— 1 



i3 



(^n ' j/n.n+g) ' Vn,n+q) 



where we define S$ = Si, as if the four NN spins were on 
a ring. In Eq. the axial and azimuthal interaction 

strengths -J f ' 9 , = D 9 ' s , -(D 9 ' s , + D 9 ' S , )/2and 

o n,n' n,n' ,zz V n,n' ,xx 1 n,n',yys' 

~ J n%[ = {Dn*n',xx ~ ^n',yy)fa as for the single-ion in- 
teraction strengths. The subscripts a,e and superscripts 
/, c correspond to our dimer notation. [5| 



C. Antisymmetric anisotropic exchange Hamiltonian 

As noted above, we write the antisymmetric 
anisotropic exchange, or Dzyaloshinskii-Moriya 
(DM), [3a, [39l Hamiltonian TL 9 njxT in the molecular 



representation, 37 



V 9 



2 6-2g 
g— 1 n— 1 



5 

n,n+q 



(12) 



We note that in these molecular coordinates, the DM 
interaction three-vectors d 9 n n+q depend explicitly upon 
the exchange bond indices n,n + q for each group g. We 
then employ the local group symmetries to relate them 
to one another. 

The rules for the directions of the d 9 n n+q were given 
by Moriya,[38j and were employed for a dimer example 
by Bencini and Gatteschi. [37[ The Moriya rules are: (1) 
d 9 n , vanishes if a center of inversion connects r n and 
tv. (2) When a mirror plane contains r n and r n >, d 9 n , 
is normal to the mirror plane. (3) When a mirror plane 
is the perpendicular bisector of r n — r„< , d? n , lies in the 
mirror plane. (4) When a two-fold rotation axis is the 
perpendicular bisector of r n — r„< , then d 9 n , is orthog- 
onal to the rotation axis. (5) When r n — r n > is an r-fold 
rotation axis with r > 2, then d g n , is parallel to r n — r n > . 
As noted above, we shall incorporate these rules in the 
molecular representation. For example, in NaX^Os, the 
lack of inversion symmetry between interacting spins has 
been shown to lead to a DM interaction. 43] 



IV. GROUP SYMMETRY INVARIANCE 
A. General considerations 

In this section, we impose the set of allowed group 
symmetry operations upon the full Hamiltonian TL 9 . 
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TABLE I: Lists of the single-ion parameter sets /xf . 

These symmetry operations are represented by the ma- 
trices 0\ for A = 1, . . . , 26 listed in Subsection A of the 
Appendix. For each g, we require H. 9 to be invariant 
under each symmetry operation 0\ for each allowed A. 

For the six g cases under study, the set {0 9 x } of group 
operations greatly reduces the number of single-ion and 
symmetric anisotropic exchange parameters. As we shall 
see, in each group <?, these reduce the single-ion and sym- 
metric anisotropic exchange interaction strength set to 

{Jj} = {Jai Je i Jf,q' Jc,q}' (13) 

for q = 1,2, which are independent of the site index n. 
That is, for each g, there are at most two single-ion, 
two NN and two NNN symmetric anisotropic exchange 
interaction strengths. In addition, for these six g cases, 
the group operations further limit the number of vector 
set parameters to 



II 9 



(14) 

i P j, (15) 

where p = q + 1 = 2,3, and we used the notation 9f = 
i = 8 9 p , etc. Some of these parameters may be 



f) 9 

"1,1' 



further restricted. In addition, however, the molecular 
single-ion and anisotropic exchange Hamiltonians contain 
both site-independent and site-dependent terms. 

For Hj-) M , we first impose the Moriya rules on each 



anisotropic exchange pair,[37J, |38| and then impose the 
required group symmetries on the six pairs. For the 
six groups under study, the group symmetries place re- 
strictions upon the d 9 n n+q , leading to the anisotropic ex- 
change parameter set 



d g 



{^zi d^i, dy ly d 9 2 , d 9 2 }- 



(16) 



For each of the six g cases, the NNN DM parameter set 
has at least one more restriction than does the NN DM 
parameter set. Some g symmetries lead to site-dependent 
signs of the components of d 9 . 



C. Imposing the group symmetries 

In Subsection A of the Appendix, we describe the ma- 
trices 0\ for A = 1, . . . , 26 representing the group sym- 
metry operations for g = C^, D^h, 0^ v , £4, D 2 d, and Td. 
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p,2 i 2°P> 3 2(p-l) 2 



TABLE II: Lists of the relevant NN (p = 2) and NNN (p = 3) 
parameter sets fif . 

9 d\ d% x d% x d 3 x2 d% 2 



C ih ,D ih d 3 z 
S 4 df 4 dlt d s y t d s j d s J 
D 2d df 2d d°? d d D J d 
T d ,C iv 
TABLE III: Lists of the DM parameter sets d 9 . 



Then, imposing Oq symmetry, we have OsS n — S n , and 



either Q\- 



tt/2 or e\ 4,1 = 0, both of which lead to 



invariance of this part of the Hamiltonian under Oq We 



therefore take the easy-axis case, 9 1 



0. Carrying 



out similar transformations on the azimuthal single-ion 
Hamiltonian leads to 



u e,n e ! 

■v Gih = di Cih 

An Vn 



Xn+1 



TT 

2' 



(23) 
(24) 



We could then choose il> < ( ih 7= 0, leaving one free angle 
parameter 4>f 4h , as listed in Table I, plus the two inter- 
action strengths j£ 4h , J^ ih . 

The symmetric anisotropic exchange Hamiltonian, 
Tiae*" 1 can be made invariant under 0\, O2, Oq in a very 
similar fashion. The operations for the other five g sym- 
metries are listed in Subsection A of the Appendix, along 
with the associated matrices. Our results for the single- 
ion, symmetric anisotropic exchange, and DM interaction 
parameters are compiled in Tables I-III. 



For each molecular group g, the allowed symmetry oper- 
ations 0\ commute with the Hamiltonian. For a partic- 
ular A, 0\r n = r„/( A ). We therefore take S n = S(r n ), 
so that OxS n = S{O x r n ) = S{r n , (x) ) = S n , (x) . 

For dh symmetry, besides the trivial identity oper- 
ation, the allowed group operations are clockwise and 
counterclockwise rotations by 7r/2 about the z axis, and 
reflections in the xy plane. [33j These operations are rep- 
resented respectively by the matrices 0i,2,6- We use this 
simple case to illustrate how the symmetries are imposed. 



D. Induced electric polarizations 



We first consider the axial part of Ti 



and set 



4 4 



n=l 



(17) 



We interpret z n *" as its vector transpose, (z ra **") T in Eq. 
([TT)) . and obtain 



OxS n 



sm9^ cos0^, 
- sin 9° ih sin cj%* h , cos d^ 4h 



(18) 



(19) 



Substituting these into the right-hand side of Eq. (TTTJ, 
setting n — > n + 1 in the left-hand side, and equating 
coefficients of S n +i ia S n +i^ for a, (3 — x,y, z leads to 



jC4h 

u a,n 


jCih jCah 

— °a,n+l ~ °a i 


(20) 




/1C4/1 r\C&h 

— °n+l — "1 ' 


(21) 




— Af ih , Z 


(22) 



As shown by Katsura et ai, 4J] the spin-orbit inter- 
actions between spins at sites n and n' can induce an 
electric polarization 



X {S n X S>i'), 



(25) 



where r n ,n' is a un it vector directed from site n to site 
n' . In our model, the thermodynamic averages of such 
polarizations vanish in the absence of DM interactions, 
but non-vanishing in-plane vector components d 9 q of the 
d 9 DM interaction parameter sets allow them to become 
finite. Tetramers with the rather low molecular group 
symmetries S4 and Z?2d have no overall center of inversion 
symmetry, and contain a complex set of DM interactions. 
Depending upon the polarizability of the attached ligand 
groups, this may lead to a combined spin- induced electric 
polarization 



1 V 

n.n' — l 



(26) 



where (. . .) represents the thermodynamic average in the 
presence of the full Hamiltonian, including the relevant 
DM interactions. 

Besides the direct DM interactions, we predict the pos- 
sibility of dual, or induced, DM interactions. Although 
d g q DM interactions between individual spin pairs are al- 
lowed in tetramers with the lowest C^h and D^h sym- 
metries studied, the group symmetry causes the dipole 
moments on opposite sides of their square geometries to 
cancel one another. Although we have not studied this 
point in detail, tetramers with these symmetries can in 
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principle be made to exhibit additional effective DM in- 
teractions by application of an electric field £/0. 44. |45 1 
Thus, in tetramers with S4, D?d, or lower symmetry, a 
multiferroic effect can occur, [45J in which both DM in- 
teractions and P s 0. More generally, multiferroic ef- 
fects arise in systems such as some dimers, trimers, and 
tetramers that generally do not have a center of inversion 
at the midpoints of the r njra '.[5( 



V. THE HAMILTONIAN IN THE MOLECULAR 
REPRESENTATION 

A. The molecular single-ion Hamiltonian 

To make contact with experiment, we use the group 
symmetries to rewrite Tt 9 f in the molecular (x, y, z) rep- 
resentation, 



2 ) 
n,y) 



(27) 



where a, = x,y,z, and we subtracted an irrelevant 
constant. TL 9 si contains the site-independent interactions 
Jf and the site-dependent interactions (— l) n J 9 and 
K^p(n), which are written in terms of the parameter sets 
fif in Subsection B of the Appendix. Most important is 
the result that for Td symmetry, 



(28) 



The first-order contributions to the eigenstate ener- 
gies from the site-dependent interactions (— f)"J| y and 
K^p(n) vanish. Hence, these interactions only contribute 
to the eigenstate energies to second and higher orders in 
the interactions Jf and Jf. For C^^D^d and S4, the ef- 
fective axial site-independent interactions Jf arise from 
a combination of the local axial and azimuthal interac- 
tions Jf and Jf . For g ^ Td, Jf can be large, even if the 
molecular structure is nearly Td- 



B. Symmetric anisotropic exchange in the molecular 
representation 

We then construct the group-invariant symmetric 
anisotropic exchange Hamiltonian in the molecular co- 
ordinates for the six g symmetries. For g = D2d,S^ 1 
there are renormalizations of the isotropic exchange in- 
teractions, modifying Tig to 



A S 2 



1 B S 



(4 



^24)1 
(29) 



J„ + SJn, 



v - v 



■SJ'g, 



where 

(30) 
(31) 

where J g , J' g are given by Eqs. ([5]) and ©, and the 6J g 

and S J g are given in in terms of the parameters sets /j,f p in 
Subsection C of the Appendix. For the three planar sym- 
metries, g = C'4h,D4h, and Ca v , 5J g = 5J g = 0. For Td 
symmetry, there are no group-satisfying azimuthal sym- 
metric exchange vectors, so Jj^ = for q = 1,2. How- 
ever, the axial vectors parallel to r n<n i satisfy all of the 
group symmetries, so that the jj d could exist, provided 

that jj\ — jji,. However, the requirement Jt a = Jj^ 
to preserve the Td symmetry of the renormalized Heisen- 
berg interactions forces jj d 2 = JJi/2. Hence, we must 

conclude that jj d q — for q = 1, 2 and 5Jr d — SJ^ d = 0. 

H%£ also leads to additional interactions SH 9 e in the 
molecular frame, 



5HL 



2 S-2q r 

EE 



Jq, z 



n,z $n-\-q,z 



+(-i) r 



T 9 

<i-! x y 



' E K9 q, a f}( n ) S n, a S n +q,fi 



(32) 



where a, [3 = x, y, z. As for the single-ion interactions in 
the molecular representation, the site-independent sym- 
metric exchange interactions J| z contribute to the eigen- 
state energies to first order, but the first-order contribu- 
tions to the eigenstate energies from the site-dependent 
interactions vanish. Both the site-independent and site- 
dependent symmetric exchange interactions are given in 
terms of the parameter sets \x\ v in Subsection C of the 
Appendix. 



C. Antisymmetric exchange Hamiltonian 

In Sees. IV and V, we already evaluated the anti- 
symmetric exchange Hamiltonians in the molecular rep- 
resentation, and the parameter sets are listed in Table 
III. These six Ti 9 DM may be combined as 



H 



1 6-2q 

dm = E E ( S » x Sn+ 9 ) ' (d 9 z (n)6 q ,iz 

q—1 n — 1 

+e^sin(mr/2) + (z X d 9 ) cos(mr/2) 

' (33) 

where the scalar d 9 z {n) and the two two- vectors d 9 q all 
vanish for g = Ci v , Td, but for the other four symmetries 
are given in Subsection C of the Appendix. 
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We note that both site-dependent and site-independent 
DM interactions give rise to second-order eigenstate en- 
ergy corrections, and can only be neglected in fits to 
experiment for tetramers with symmetries very close to 
TdjC-to, or higher. In Subsection D of the Appendix, we 
give H 3 DM for the lower symmetry tetramers. 



D. Biquadratic and three-center quartic isotropic 
exchange interactions 



the six g symmetries may be written as 

n = XX, 

9=1 



n -t,l 



K2 



(36) 



9 _ 



4 4 



'^t,l XZ XZ C^ n ' ^n+l){S n > ■ Sn'+l)) (37) 

n=l „/=i 

odd cvcn 
4 2 

~Jt,2 2J a2 ' Sn+l) (S n > ■ S n > +2 ). (38) 



n— 1 n' — 1 



In the previous subsections, we listed the quadratic 
single-ion and anisotropic exchange interactions for the 
six high-symmetry tetramer groups under study. How- 
ever, in the lower-symmetry AFM tetramers {M4M012}, 
with C\ v symmetry, [46[ and M4 and C04 [2 x 2] grids (or 
rhombuses), with approximate symmetry. [471 l48l 1-491] 
fits to magnetization data were facilitated by the inclu- 



sion of biquadratic interactions. [46J, |47j, |48j In the for- 
mer case, the powder fits assumed field-dependent inter- 
action parameters, but the single-ion interactions were 
assumed to have symmetry, which vanish to first or- 
der in their strength, and the anisotropic exchange inter- 
actions were neglected. In the latter case, the au- 
thors neglected the NN DM interactions given in the Ap- 
pendix. Subsequently, Kostyuchenko showed that three- 
center isotropic quartic interactions should be compa- 
rable in magnitude to the biquadratic interactions, and 
provided a fit to the midpoints of the level-crossing mag- 
netization behavior on {M4M012}, without making the 
assumption of st rong field dependence to the Heisen- 
berg interactions. [50j Here we provide preliminary fits 
to the AFM {M4M012} magnetization data, extending 
the treatment of Kostyuchenko to include the first or- 
der anisotropy interactions, which can fit the widths of 
the transitions, as well as the midpoints. More com- 
plete fits to those experiments and to ex per iments on the 
grid SMM's will be presented elsewhere. 5l| Such fits are 
greatly aided by an analytic treatment of biquadratic and 
three-center isotropic quartic exchange. 

For tetramers with the six g symmetries under study, 
the biquadratic interactions may be written as 



Ti 9 



6-2q 
n=l 



(34) 
(35) 



J b 2, but otherwise J 9 2 ¥" 



For g = T d , we take 1 

For the six g symmetries, TL^ is invariant under all 
of the appropriate symmetries. 

The three-center quartic interactions for systems with 



VI. EIGENSTATES OF THE FULL 
HAMILTONIAN 

A. Induction representation 



We assume a molecular Hamiltonian 

n 9 = HV + ni 



XU9 1 n-/9 

urL ae ' rL DM 



(39) 



To take proper account of B in 7io' r , we construct our 
SMM eigenstates in the induction representation by 









-( 








— sin <p sin V cos < 
cos 4> sin 9 sin < 
cos6> 



V 

z' 



(40) 



so that B = Bz' . A subsequent arbitrary rotation about 
z' does not affect the eigenstates. [j| We then set h = 1 
and write 



s 2 IV^' S24 ) 




(41) 




= S13(S13 + 1M^ 24 ), 


(42) 




= S2 4(S24 + 1M^ 24 >, 


(43) 






(44) 


'S'o-IV'a.m' 824 ) 


= A™\^i), 


(45) 


A™ 


= y{s — m)(s + m + 1), 


(46) 



where S& = Sx + i&Sy with a = ±. For brevity, we define 



v = {s,m,si 3 ,s 2 4,si}, 
v = {s, sis, S24, si}, 



(47) 
(48) 



where V excludes m, and write \v) = |V4.m' S24 )- 
and 7io' r are both diagonal in this representation, but 
TCsijS'Hae, and l are not, and we therefore assume 
their interaction strengths to be small, relative to J g 
and J' g . From Eqs. @T]) to gD, {v>\H% r + H 9 hfi \v) = 
qS u i_ u , where 



E 9 



Jo 



s(s + 1) - "fBm + 5EjL Q 



(49) 
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, 2 

^1,0 = - 2 X! [(4 ~ jg) S n,n+2(Sn,n+2 + 1) 



J9 2 

— ^-[-2si(si + 1) + s n<n+2 (s. 



n,n+2 



i)] 2 



i 9 

J t,2 i 



-I — 2 J -[-2si(si + 1) + S n ,n+2(Sn,n+2 + 1)] 



X fs(s + 1) - ^ s «',n'+2(s«%n'+2 + 1' 



n'=l 



(50) 



where the J s and are given by Eqs. (JSJ), ©, and 
(|129[) - (|132[) . Since 7ig' r and 2 are invariant under all 
rotations, E 9 is independent of 8, <fi. 
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B. First-order eigenstates 

In the induction representation, we write "H%i + 67i% e + 
H 9 DM as H' s 9 + 7i[i 9 e + H'p M ■ W£ x is a scalar independent 
of the direction of B. We then make a standard per- 
turbation expansion for the seven remaining microscopic 
anisotropy energies {J?} = {J 9 , J 9 , J 9 f q , J 9 q , J^J for 

g = 1,2 small relative to |J 9 |, |Jg|.f$] To do so, it is nec- 
essary to evaluate the single-ion matrix elements analyt- 
ically, as they contain much of the interesting physics. 
Compact expressions for these matrix elements for gen- 
eral (si, S2, S3, S4) are given in Subsection E of the Ap- 
pendix. 

At arbitrary B angles (8, </>), the first order corrections 



E 9 



TL ' 9 



rL DM 



V 9 



n? 



to the 



eigenstate energies for g = C^h, D^h, C^ v , S4, Did, and 
Td symmetries are 



E 9 



j?- v 



[m 2 - s(s + 1)] - 5 J 9 ' " 



[3TO2 -f + 1)] jrcos^, (51) 
' (52) 



\ J IM- 



(53) 



where analytic expressions for the oi , &i , ci, #77 and 7^- 
for general 77 are given in Subsection F of the Appendix, 
along with Tables IV- VII and VIII-XI of their simple 
analytic forms for the lowest four eigenstate manifolds 
of FM and AFM tetramers, respectively. We note that 
all of these interaction coefficients are invariant under 
S13 <-* S24, as expected. The DM and all site-dependent 
interactions vanish in this first-order perturbation. Sec- 
ond order corrections to the eigenstate energies will be 
presented elsewhere. 51 



FIG. 4: Plots of the magnetization M/7 versus yB/\J g \ of an 
si = 1/2 tetramer at T — with c = a, g — Si, d g z = d\ y = 

d L = d 2y = rf . d lx = °. d /\Jg\ = 0-1, and B||(lll). The solid 
and dashed curves are for the Type I (J' g /J g = 1.5) and Type 
II (J'g/Jg = 0.5) tetramers, respectively. 



For all six g symmetries, E 9 1 has a form analogous 
to that of the equal-spin dimer in the absence of az- 
imuthal single-ion and symmetric anisotropic exchange 
interactions. |5| For these high-symmetry tetramers, to 
first order in the anisotropy interactions, the azimuthal 
single-ion and anisotropic exchange interactions merely 
renormalize the respective effective site-independent ax- 
ial interactions. Thus, to first order, we only have two 
effective isotropic exchange, two biquadratic isotropic ex- 
change, and three effective anisotropy interactions, J ff , 
J' g , J 9 17 J 9 2 , J 9 , J 9 Z , and J| z , which are fixed for a 
particular SMM. Nevertheless, the first-order eigenstate 
energies E 9 = E 9 + E 9 v l , given by Eqs. (gSJ and lj5T]). 
contain these seven effective interaction strengths in ways 
that depend strongly upon the quantum number set v 
and upon 8. These different i/, 8 dependencies can be 
employed to provide definitive measures of at least some 
of the seven 77-independent effective isotropic exchange 
and anisotropy interactions. 



C. Type I and Type II tetramers 

There are at least two types of FM and AFM 
tetramers. To the extent that single-ion, symmetric 
anisotropic exchange, and biquadratic exchange interac- 
tions are small relative to the Heisenberg interactions, 
there are just two types of tetramers. The criterion is 



simply based upon J' — J g 



in SE 9 



which we assume 



to be larger in magnitude than all anisotropy and bi- 
quadratic interaction strengths. Type I tetramers have 
J' — J„ > 0, which can occur for either sign of </„, pro- 
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FIG. 5: Plots of the magnitude of the spin-derived polariza- 
tion |P S | in arbitrary units versus jB/\J g \ at T = of an 
si = 1/2 tetramer with c = a, g — S4, d 9 z = df y — d% x = 
d 9 2y = d, d{ x = 0, d/\J g \ = 0.1, and B||(lll). The solid and 
dashed curves are for the Type I (J' g /Jg = 1.5) and Type II 
(J'g/Jg = 0.5) tetramers, respectively. 
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FIG. 7: Plots of the magnetization M/7 versus ^B/\J g \ of an 
si = 1 tetramer at T = with c = a, g — S4, d% 



d, d 9 lx = 0, d/\J g \ = 0.01, and B||(lll). The 



solid and dashed curves are for the Type I {J' g /J g 
Type II (J' g /J g = 0.5) tetramers, respectively. 
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FIG. 6: Plots of the magnitude of the spin-derived polar- 
ization |P S | in arbitrary units versus 9/ir at T = of an 
si = 1/2 tetramer with c — a, g — S4, J'g/Jg = 1, 
d% = d\ y = d\ x = d% = d, d\ x = 0, d/J g = 0.1, and <f> = n/4. 

The solid and dashed curves are for r yB/\J g \ = 1.5, 1, respec- 
tively. 



vided g ^=Td- For Type I, the lowest energy state in each 
s manifold occurs for the maximum values S13, S24 = 2s 1. 
Thus, at low T, Type I tetramers behave as pairs of spin 
2si dimcrs. Type II tetramers with J' g — J g < are 
frustrated, with the lowest energy state in each s mani- 
fold occurring for the minimal si3,S24 values. For even 
s, these minima occur for si3,s 2 4 = s/2, but for odd s, 
the energy minimum is doubly degenerate, occurring at 
Si3,s 2 4 = (s ± l)/2, (s=F l)/2- Hence, explicit formu- 



FIG. 8: Plots of the magnitude of the spin-derived polariza- 
tion |P S | in arbitrary units versus jB/\J g \ at T — of an 
si = 1 tetramer with c = a, g — S4, d% — df y — d\ x = 
d\ y = d, d\ x = 0, d/\J g \ = 0.01, and B||(lll). The solid and 
dashed curves are for the Type I {J' g /J g = 1.5) and Type II 
[J'g/Jg = 0.5) tetramers, respectively. 



las for first-order cigcnstate energy parameters a^, b£, 
ci, and Bjt with arbitrary s in the three special cases of 
(si 3 ,s 24 ) = (2si,2si), (si 3 ,s 24 ) = (s/2, s/2) for even s, 
and (S13, S24) = [(s± l)/2, (s =p l)/2)] for odd s are given 
in Subsections G and H of the Appendix. For sufficiently 
strong J' g — J g , these formulas can apply to the lowest 
energy eigenstate in each s manifold. When J' g — J g is 
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0.50 



FIG. 9: Plots of the magnitude of the spin-derived polariza- 
tion |P 8 | in arbitrary units versus 9/tt at T = of an si = 1 
tetramer with c = a, g = £4, J'g/Jg = 1, d% = d\ = d\ x = 

a 2y 



d, d{ x = 0, d/J g = 0.01, and </> = tt/4. The solid and 



dashed curves are for ■yB/\J g \ = 1.5, 1, respectively. 



small relative to the other interactions, the situation be- 
comes more complicated, as the lowest energy eigenstate 
for a particular B can depend upon more of the effec- 
tive interaction values. Tables IV- XI in the Appendix 
are sufficient for full analyses of such cases for s\ < 3/2, 
but such studies will be made subsequently. 51 1 



D. DM interactions and spin-induced polarizations 

Although the first order correction to the energy aris- 
ing from the DM interactions vanishes, for s± = 1/2, 1, it 
is not too difficult to diagonalize the Hamiltonian matrix 
exactly, and hence to take precise account of the effects of 
the DM interactions. To focus on the effects of DM inter- 
actions, in Figs. 4-9, we omit the single-ion, symmetric 
anisotropic exchange interactions, and biquadratic inter- 
actions, keeping only the AFM Heisenberg, Zeeman, and 
weak DM interactions. In Fig. 4, we plotted M/7 ver- 
sus -fB/\J g \ for si = 1/2 AFM tetramers with £4 sym- 
metry with B||(lll), T = 0, and d/\J g \ = 0.1, where 
d 9 = d\ y = d\ x = d 9 2y = d, d\ x — 0, in the limit 



d 9 

a 2y 



c = a. We note that for J g < and J' g /J g = 0.5 the 
magnetization exhibits sharp steps close to the integral 
values s = 1, 2 of jB/\J g \, characteristic of Type I AFM 
tetramers. For J g < and J' g /J g — 1-5, the two steps are 
shifted to higher B values, and the first step is broadened. 
This is Type II AFM tetramer behavior for the DM inter- 
action. In Fig. 5, the corresponding curves for the spin- 
induced polarization \P S \ are shown. In both cases, there 
is a sharp peak at the inflection point of the first magne- 
tization step, at which the total spin value s changes from 



to 1. In addition, there is a discontinuity in slope at 
the positions of the second magnetization steps, at which 
P s begins its rapid decrease to zero, which it reaches at 
the B value for which M just reaches saturation. Note 
that P s — > at large B, since the large B aligns the 
spins at all sites, causing their vector products to van- 
ish. Some results for the intermediate case J' g /J g = 1 
are shown in Fig. 6. In this figure, we plotted \P S \ in 
arbitrary units versus 9/tt for s\ = 1/2 AFM tetramers 
with S4, symmetry at cj> — 7r/4 with d/\J g \ = 0.1 with 
7-B/l Jg I = 1-5,1, respectively. The curves are symmet- 
ric about 9 = 7r/2. We note that the curves exhibit 
broad maxima at 9/tt w 0.3 for both B values. These 
broad curves reflect the strong frustration of the spins 
with J' g = J g , for which no preferred s%3, S24 values exist. 

In Figs. 7-9, we plotted analogous curves for si = 1 



tetramers, except that d/\J g 



0.01. Curves with 



d/\J g \ = 0.1 for si = 1 exhibit steps or peaks which 
are much more broadened than the corresponding ones 
for s% = 1/2 pictured in Figs. 4-7. We note that in Fig. 
7, the dashed curve for the Type II case J' g /J g — 1.5 has 
a wider steps at M/7 = 0, 2 than at 1,3. This low value of 
the DM strength d/\J g \ = 0.01 leads to very sharp peaks 
in \P S \ for both Type I and II tetramers, as shown in 
Fig. 8. However, the three peak tails are much broader 
for Type II tetramers with J' / J g = 1.5 than for Type I 
tetramers with J' g /J g = 0.5. In each case, the peak po- 
sitions correspond to the first three magnetization step 
jB/\J g \ values, and the polarization also vanishes at the 
jB/\J g \ value at which the fourth magnetization step is 
completed. Finally, in Fig. 9 we plotted \P S \ in arbitary 
units versus 9/n with the same parameters as in Fig. 6, 



except that d/\J g \ 
has a peak at 9/w 



: 0.01. The curve with jB/\J g \ = l 
0.3, as for the corresponding curve 
with s% = 1/2, but the solid curve for ~fB/\J g \ = 1.5 
has a flat region for 9/ir between 0.3 and 0.45, and both 
curves become vanishingly small at 9 — n/2, which differs 
strongly from the behavior shown in Fig. 6 for si = 1/2. 
This difference suggests an interesting parity effect at 
9 = tt/2, with P s (9) oc (0 - tt/2) 2si+1 . 



E. First-order AFM level crossing inductions 

For AFM tetramers, J g < 0. There will be 2s 1 + 1 level 
crossings, as exhibited by the magnetization steps in Figs. 
4 and 7 for si = 1/2,1, respectively, provided that the 
lowest energy state in each s manifold does not exhibit 
level repulsion. In order to specify the level-crossing in- 
ductions, we first write E§ m (si3, S24, si) = E 9 + E 9 1 . 
We then note that the si3,S24 values involved in level 
crossings are those corresponding to the lowest energies 
for a particular s value. These are different for Type I and 
II tetramers. For Type I tetramers, the level crossing- 
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inductions occur for 



si = 1/2 first-order AFM level crossings 



£f )S (2si,2 Sl ,si) = £f_ M _ 1 (2 Sl ,2si, Sl ). (54) 

and for Type II tetramers, they occur for 

El s (s/2,s/2, Sl ) = E 9 s _ hs _ 1 [( S -l±l)/2, 

(a - 1 =F 1)A «i] (55) 



for even s, and 



E 9 

S.S 



( S ±l)/2,( ST l)/2, Sl ] = El 



(s-l)Asi 



l)/2, 
(56) 



for odd s. In Subsections G and H of the Appendix, we 
presented the formulas for the level-crossing induction 
parameters for both Types I and II tetramers. 

For Type I tetramers, the first-order level-crossing in- 
ductions B 9 ^ 1 ^ obtained from Eq. (f5"4")l have the re- 
markably simple form, 



For the simplest case S\ = 1/2, as in AFM CU4 
tetramers, the single-ion interaction Jf does not con- 
tribute to the level-crossing inductions, as for the dimcr 
of equal s\ — 1/2 spins. [5j Using the results given in 
Subsections G and H of the Appendix, the expressions 
for the 7-8^2 ^ (9) functions are particularly simple. For 
,si = 1/2 effective-dimer Type I tetramers, J' g — J g > 0, 



+ — *^ 2 ) 

+~J e yi-3cos 2 0), 
-2 7 4- J 9 - T 9 

-ij o 9 ff (l-3cos 2 0), 



(60) 



(61) 



with 



= 1/2. For 

frustrated Type II tetramers with s% = 1/2, J g — J ' > 0, 



where JfL is given by Eq. 



lB 9 - l f>(9) = -sJ g 



b+ 

J9_ J9 A 

J z 2 ' ' 



3 ( J ° s 



(2si - 1) 



(4si - 1) 



Jf (1-3 cos 2 



T 9 



J 9 S! J 9 

J i,z , 



2.z 



4si - 1 ' 



(57) 
(58) 



where 6 



^/ 1+ ( s )' c i = c Ti ( s )' ^ — (^/(s), and 



e = ej 1 (s) are given in Subsection G of the Appendix. 
We note that this is independent of J' g and J 9 2 , and that 
the NN and NNN symmetric anisotropic exchange inter- 
actions combine to yield the universal Type-I level cross- 
ing form J 9 S . Furthermore, the ^-dependencies of the 
single-ion and symmetric anisotropic exchange contribu- 
tions have the same s-dependencies for fixed s\. However 
the ^-independent contributions from Jf, J 9 S , and J 9 1 
depend separately upon s for fixed s%. 

For the Type II AFM tetramer level-crossing in- 
ductions, the contributions from the near-neighbor 
anisotropic exchange interaction Jf has a rather simple 
form. As shown in Subsection H of the Appendix, these 
contributions 7^;' c(1) = J 9 J ltZ {s,9) to 7 Sf; k s (1) (9) are 
independent of s\, where 



^(l + (2 S -l)cos 2 ( 
3( ^ TT (l + (2.s-3)cos 2 0), 



s odd 



(59) 



s even. 



Note in particular that for s = 1, /i, z (l,0) = 0. How- 
ever, the single-ion and NNN symmetric anisotropic ex- 
change contributions to the level crossing inductions de- 
pend upon both s and s± in different ways. 



jB^(9) 



+ ^ 2 yi + cos 2 9), 
- - T - V 4- -( T 9 - J 9 \ 



(62) 



"2(^,2 ^2) 



■■9). (63) 



Even in this simplest of all tetramer cases, there is still 
a qualitative difference between the level crossing induc- 
tions of Type I and Type II AFM s\ = 1/2 tetramers. For 
Type I, there is only one effective anisotropic exchange in- 
teraction, J 9 S = (Jf Z +J 9 z )/2 that affects the level cross- 
ing. However, for Type II s\ = 1/2 tetramers, the level 
crossing is different for the NN and NNN anisotropic ex- 
change interactions. The only effect of the group symme- 
try is to provide restrictions upon the values of the inter- 
actions. These expressions also show that Type II AFM 
tetramers have a more complex level-crossing induction 
variation than do Type I AFM tetramers, as the first- 
order Type I level-crossing behavior is fully described by 
three parameters, whereas the first-order Type II level- 
crossing behavior depends upon four independent param- 
eters. On the other hand, for this special s\ — 1/2 exam- 
ple, the 9 dependencies of the first and second jB 9 '^^ 
are opposite in sign for Type I, but can have the same 
sign for Type II. In Fig. 10, we illustrate these si = 1/2 
behaviors for the Type I with J 9 s /J g = 0.2 and for Type 
II with 3 9 2 JJg = 0.2 and J g - J' g = 0.5| J g \. 

In the special case of T4 symmetry, we have Jx d — 

Jbt, and 



J b,2 



Jb, Jt.l 



• J t,2 
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FIG. 10: Plots of the si = 1/2 first-order level-crossing 
7Bji°' g (0) /\Jg\, where 8 is the angle between B and z, with 
< q = 0, J c 9 ff /J 9 = 0.2, J' g - J 9 > (solid, Type I) and 
JlJJg = 0.2, J 9 - J' g = 0.5 | (dashed, Type II). 

jji, = 0. Since for s\ — 1/2, Ti. 9 b and are diagonal, 
the Jb and J* dependencies of the eigenstate energies are 
exact. Hence, the only difference between Type I and 
Type II si = 1/2 tetramers with Td symmetry is deter- 
mined by the sign of J& — Jf, as is evident by comparing 
Eqs. (fuT)f - (163p . In any event, there is no 9 dependence 
to the first-order level-crossing inductions for s% = 1/2 
AFM tetramers with Td symmetry. 

si > 1 first-order AFM level crossings 

For si > 1/2, the AFM level-crossings become much 
more complex than for the Si = 1/2 case, as single-ion 
anisotropics Jf are allowed, and the biquadratic inter- 
actions J^ affect the various level-crossings differently. 
We first consider the simplest si > 1/2 case, s% = 1, 
appropriate for AFM N14 tetramers. Exact expressions 
for the s = 1,2,3,4 first-order level-crossing inductions 
Bf' l ^ 1 \0) for Type I and II Si = 1 tetramers are given 
in Subsection F of the Appendix. In Figs. 11 and 12, we 
plotted the 0-dependence of the first-order level crossing 
induction -yBf^ {9)/\J g \ for g = C Ahl D ih , C iv , S 4 , and 
D 2 d, for two AFM Type I examples and for three Type 
II examples with J g — J' g = 0.5| J g \, respectively. In each 
curve, we allow only one of the anisotropy interactions 
(or effective interactions) to be non- vanishing. In Fig. 
11, the solid and dashed curves are for J^/J g = 0.2 and 
J^s/Jg — 0.2, respectively, where J% s = Jf z /2 + J| z /3 
for si — 1. We note from Eq. (f57|) and from Fig. 11 
that for Type I, the single-ion and symmetric exchange 
anisotropies lead to opposite ^-dependencies, both hav- 
ing a change in sign just before the second level crossing, 
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FIG. 11: Plots of the first-order level-crossing B^ (1) (8) /\j g \ 
for Type I, J g — J' g < 0, J^ q = 0, and si = 1. Solid curves: 
Jg/Jg = 0.2. Dashed curves: J 3 cS /J g = 0.2. 

and the dependence of 'yBf^ X ' upon 1 decreases with 
increasing s. 

In Fig. 12, we illustrate the si — 1 Type II level cross- 
ings, setting J g — J' g = 0.5| J g \. The solid curves are for 
Ji/Jg = 0.2 for g = Cy*, £>4fc, C iv , S4, and D 2d , as in Fig. 
11. The dashed and dotted curves are for Jf Z /J g = 0.4 
and Jf Z /J g = 0.4, respectively. For each curve, the Type 
II isotropic exchange parameters lead to a larger gap be- 
tween the s — 2 and s = 3 level crossings. The sign of the 
^-dependencies of the single-ion (solid) curves changes 
between s = 2 and s = 3. The effects of the NN symmet- 
ric anisotropic exchange interactions vanish for s = 1 , but 
increase in magnitude with increasing s for s = 2,3,4. 
The sign of the ^-dependence of the level crossing due 
to the NNN symmetric anisotropic exchange interactions 
does not change, but its magnitude increases monotoni- 
cally. Thus, Type II AFM si = 1 tetramers have a richer 
set of first-order level-crossing behaviors than do Type I 
AFM s± — 1 tetramers, and the unpictured biquadratic 
interactions increase this richness. 

In Figs. 13 and 14, the analogous Type I and Type 
II first-order AFM level crossing inductions are plotted 
versus 9 for si = 3/2 equal-spin tetramers, such as Cr±. 
The notation is the same as in Figs. 11 and 12. For 
Type I Si = 3/2 AFM tetramers, the single- ion and sym- 
metric anisotropic exchange interactions lead to different 
^-dependencies of the level-crossing inductions, each with 
a change in sign in the 9 dependence at about the second 
level crossings, as seen in Fig. 13. For Type II si = 3/2 
AFM tetramers, the sign changes appear between the 
second and third level crossings, as shown in Fig. 14. 
Although not pictured, the contributions for si = 3/2 to 
the level-crossing inductions from the can be easily 
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FIG. 12: Plots of the first-order level-crossing 
~fBlf° {V >{e)/\J g \ for Type II with g = C Ah , D Ah ,C iv , Si, D 2d , 
J g — J' g — 0.5\J g \, J^ q = 0, and si = 1. Solid curves: 
Ji/Jg = 0.2. Dashed curves: Jf Z /J g = 0.4. Dotted curves: 

JlJjg = 0.4. 
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jB*$?(6)/\j„\ for Type I, J g - J' g < 0, J b s q = 
and si = 3/2. Solid curves: Jf/J = 0.1. Dashed curves: 



Plots of the first-order level-crossing 
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= 3/2. Solid curves: Jf/J 9 
0.1. 



calculated from Eq. (|50|) and the expressions for ^//^(s) 
and d 3 /j 2 e (s) in Subsection H of the Appendix. They con- 
tribute to a more complex level-crossing pattern than for 
a x = 1. 

For g — Td with s\ > 1/2, we still have J| 2 = 
Jf = 0, so that the first-order level-crossing inductions 
75j,' s lc ' 1 ' {0)/\ Jr d | vary from integral values only by 6- 
independent constants due to the biquadratic interaction 



VII. THE SELF-CONSISTENT HARTREE 
APPROXIMATION 

A. Partition function and thermodynamics 

The self-consistent Hartree approximation, or strong 
exchange limit, [13] provides accurate results for the B 
dependence of the specific heat and magnetization at 
low kBT/\Jg\ and jB/\J g \ not too small, [H where ks 
is Boltzmann's constant. In this approximation, Ef, = 
E g vS) + E g vl is given by Eqs. ((49|) and (|51~j) . respectively. 
We shall present the self-consistent Hartree approxima- 
tion of four measurable quantities in the induction rep- 
resentation. We first define the trace valid for our eiqen- 
state representation, 

2si S S13 + S24 

T ^-E= E E E > ( 64 ) 

V Sl3,S24=0 771=— S s=|«13— S24 

The partition function in the self-consistent Hartree ap- 
proximation may then be written 



Z<?> = Tr u e 



:->E'> 



(65) 



where (3 — l/(fcsT). In this compact notation, the self- 
consistent Hartree magnetization Mg X \B,T) and spe- 
cific heat Cqy(B, T) are given by 



M^(B,T) = 7 Tr„( 



(66) 
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We note that there are strong differences between the 
low-T behavior of FM and AFM tetramers. We assume 
\Jg\ > \J'g~ Jg\- For FM tetramers with J g > 0, the low- 
T thermodynamic behavior is dominated by the s = 4si, 
m = — 4si state, leading to 

M«(B,T) T ~ Q jBB^iPlB), (68) 

where Bs(x) is the Brillouin function. The universal- 
ity of this function renders thermodynamic studies use- 
less for the determination of the microscopic parameters. 
For AFM tetramers with J g < 0, however, there will 
be interesting level-crossing effects, which can be em- 
ployed to measure the microscopic interaction parame- 
ters, as discussed in detail in Sec. V. As for dimers, 
C V (B,T) for AFM tetramers at sufficiently low T ex- 
hibit 2si central minima at the level-crossing inductions 
B l'us(9) ~ B 9 s^s {1) {9) that vanish as T -> 0, equally sur- 
rounded by peaks of equal height. 0] As for the magne- 
tization, Cv(B,T) for FM tetramers at low T reduces 
to that of a monomer with spin 4si, yielding a rather 
uninteresting Schottky anomaly. 

B. Electron paramagnetic resonance 

However, the microscopic nature of FM tetramers can 
be better probed either by EPR or INS techniques. The 
self-consistent Hartree EPR absorption S$X-a a(B , w) for 
clockwise (a = 1) or counterclockwise (er = —1) circu- 
larly polarized oscillatory fields normal to B is 

3y 9 ' (1) - -^-Tr Tr ,e~^\M >\ 2 

x [6(Ei - El, + w) - 5{El, El + w )] , 

(69) 

where M v y = A^ m S m ^ rn+a S s > iS 6 s ' i3ySi3 S s > 24 ^ 24 and TV = 
Y] ,. The strong resonant inductions appear at 

7 BS.« = ± w + ^±^(l-3cos^)Jr, 

(70) 

where Jf 17 is given by Eq. flS2J| We note that if- 77 
contains the three effective microscopic interactions, Jf, 
Jf z , and Jf z , multiplied by the constants ai, — c—, and 
— a^/2, respectively. In Tables IV and VI of Subsection 
F of the Appendix, the values of these parameters for 
the FM ground state and the first three excited state 



manifolds for arbitrary si are given. We note that EPR 
measurements are insensitive to the Heisenberg and bi- 
quadratic exchange interactions, which preserve m. 

For either FM or AFM si = 1/2 tetramers, EPR mea- 
surements can only probe the two microscopic symmet- 
ric anisotropic exchange interaction parameters Jf z and 
Jf , and measurements of the two s — 1 excited states 
are sufficient to determine them. For FM tetramers with 
si > 1, it is a bit more difficult. From Tables IV and VI 
and from Eq. (|219[> in Subsection G of the Appendix, it 
is easily seen that the (s, sia, S24) = (s, 2si, 2si) states all 
provide measurements of the same combination of these 
three microscopic interactions. Hence, for FM tetramers 
with si > 1/2, measurements of the ground s = 4si and 
the first excited state manifold with s = 4si — 1 are in- 
sufficient to completely determine the three microscopic 
interactions. In order to stay within a single s value for 
FM tetramers, one would need to study the second (or 
higher) excited state manifold with s = 4si — 2 (or lower), 
in order to obtain sufficient information to determine 
the three microscopic interaction strengths. For AFM 
tetramers with s\ > 1, EPR transitions in the ground 
state are not allowed, but measurements of the first ex- 
cited s = 1 state manifold would suffice to determine Jf 
and the J qz , as seen from the formulas in Tables VIII 
and X in the Appendix. 

C. Inelastic neutron scattering 

The Hartree INS cross-section Sg 1 ^ (B, q, 10) is 

x5(<j + E° -El,\ (71) 

where a, (3 = X,y,z, g £ = sin iig cos % = 
sin 6b,g sin <f>b,q , an d qg — cos0b l9 , 6b „ and <pb,q describe 
the relative orientations of B and qr,[5j the r n are given 
by Eq. (P), and the (v'\S n ^\v) are given by Eqs. (|164[) 
and (|165p in Subsection E of the Appendix. The scalar 
q ■ (r n — r„<) is invariant under the rotation, Eq. (|40|) . 
After some algebra, we rewrite Sg 1 " 1 (q, to) as 

v' 

x (sin 2 e h , q L vy {q) + ^^^(?)) , 

(72) 

where the Hartree functions L u y(q) and M u y(q) are 
given in Subsection I of the Appendix. They are inde- 
pendent of B. Since E® is well-behaved as B — > 0, Eq. 
(TT2")) is accurate for all B. 
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As for the dimer,Q additional EPR and INS transi- 
tions with amplitudes higher order in the anisotropy pa- 
rameters Jf, Jf z , and J| z relative to J g are obtained 
in the extended Hartree approximation, but will be pre- 
sented elsewhere for brevity. [BH 



VIII. DISCUSSION 

The quadratic phenomenological total spin anisotropy 
model widely used in fitting experimental data on SMM's 



H,, 



A-VS 2 -£(S 2 x -S 2 y ), 



(73) 



where A represents the isotropic total spin interactions, 
and T> and £ are measures of the axial and azimuthal 
total spin anisotropy, respectively. Q] Often, additional 
quartic terms are added. jH, The anisotropy is de- 
fined relative to the total spin principal axes, which for 
equal spin, high-symmetry systems are the molecular axis 
vectors. It is easy to evaluate = {v\H p \v) in the in- 
duction representation. One obtains Eqs. (|49|) and (f5Tj) . 
provided that 



A 
Ao 
6A 
V 
£ 



Ao + SA, (74) 

-J g s{s + l)/2~ 1 Bm, (75) 

SE* fi - SJf, (76) 

ir, (77) 

0, (78) 



where 6E% fi , Jf v and SJ%' V are given by Eqs. (50]), (5J) 
and (|53|) . respectively, which contain the constants a^, 
b^, ci, and B-. Precise formulas for arbitrary 77 ap- 
pear in Subsection F of the Appendix, along with Tables 
IV- XI of their values for arbitrary s\ in the ground and 
lowest three excited state manifolds for FM and AFM 
tetramers, respectively. Usually, one assumes the strong 
exchange limit, so that the isotropic Ao is sufficiently 
large that it remains constant for B = 0, and can be ne- 
glected. A non-vanishing £ would lead to a term in Eq. 
(f5Tj) proportional to sin 2 8 cos(20), as for the dimer,[5| 
which does not arise in the first-order calculation for the 
high-spin tetramers under consideration, based upon the 
microscopic parameters alone. Hence, the T> term in Ti. p 
alone describes the 6 and m dependencies of E a v + E 9 V 1 
correctly, provided that the quantum numbers S13 and 
S24 remain constant. 

More important, the additional constant term SA has 
generally been neglected. Even to zeroth order, the sign 
of J g — J g in SE^L distinguishes between Type I and 
Type II tetramers, which distinction is absent in the 
phenomenological model. Moreover the different first- 
order dependencies of A and T> upon the V are impor- 
tant in determining the level-crossing inductions for AFM 



tetramers, each of which involves two values of s and m, 
and hence different S13 and S24 values, as well. The zero- 
field energy spectrum is thus more complicated than that 
given by the usual phenomenological model, which could 
lead to substantially different fits to experiment. 

For simplicity, the only higher order interactions we 
have considered are the isotropic NN and NNN bi- 
quadratic exchange interactions. These isotropic interac- 
tions are rotationally invariant, so they are independent 
of 6 in the induction representation. Hence, they only 
contribute to SA. Thus they modify the positions but 
not the ^-dependencies of the AFM level-crossing induc- 
tions, and do not modify any EPR transitions. 

In the ground state of FM tetramers, V is restricted to 
the single set of values, (s, S13, S24) = (4si, 2s\, 2si). In 
this high-spin case, TL p can provide a correct phenomenol- 
ogy of the ground state energy. However, if applied to the 
two low- lying excited states with s — 4si — 1, for instance, 
one would infer two different A and T> values from those 
obtained in the ground state. However, as noted above, 
since all states with (si3,S24) = (2si,2si) contain the 
same combination of and c^, in order to exploit the 
v dependence of T> to obtain an unambiguous EPR mea- 
surement of the three microscopic parameters Jf, Jf z , 
and J| 2 , for si > 1/2, one needs to examine higher state 
manifolds, such as the manifold with s = 4si — 2. For 
AFM tetramers, H p also correctly provides a vanishing 
first-order correction to the s — manifold of states with 
S13 = S24 = 0, 1, . . . , 2s\. However, 7i p also has problems 
describing the first excited manifold of AFM states with 
s = 1, because it leads to the choice of 8 A independent 
of the quantum numbers v, which is unphysical except 
for tetramers with Td symmetry. Hence, the phenomeno- 
logical model works best in describing only a single state 
with fixed (s, S13, S24). This is more restrictive than the 
usual assumption of its applicability to all states with 
fixed s.fHillUlHlH! 

We note that the FM Cu 4 tetramer Cu 4 OCl 6 (TPPO) 6 
was claimed to have Td symmetry and an s = 2 
ground state. 12, [IH, 14 1 It is noteworthy that those 



authors thought that anisotropic exchange interactions 
might be responsible for their observed zero-field energy 
splittings. [l2L [l5| Since tetramers with Td symmetry do 
not have either symmetric or antisymmetric anisotropic 
exchange interactions, another explanation must be con- 
sidered. From Tables IV- VII in the Appendix, it is ev- 
ident that the FM ground state is non-degenerate for 
all si, even for those SMM's with lower symmetries 
allowing anisotropic exchange interactions. It there- 
fore appears that the sample may not have been single- 
phase, 12J,|l3, [l4 1 as in a nominally S4 M4 tetramer, [22j 
allowing for an apparent ground state splitting. 

We note that for the FM Fe 4 SMM, Fe4(thme) 2 (dpm) 6 , 
where Hsthme is 1,1, 1-tris (hydroxymethyl)ethane and 
Hdpm is dipivaloylmethane, 52|, l53l . |54| | the high D3 sym- 
metry also precludes the £ term in 7i p . Nevertheless, 
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in fits to INS data, it was assumed that £ ^ 0, H2, 55 1 
in order to obtain the appropriate anticrossing gaps, so 
that either the powdered sample did not have pure D 3 
symmetry, or the phenomenological model they used, Eq. 
([75)1 plus two quartic terms obeying D 3 symmetry, was 
not appropriate. Since single-ion interactions appeared 
to be important .(5^1 the total spin might not have been 
a well-defined quantum number, as in at least one Fe2 
dimer and in Feg.0, d, S, Q However, we note that it 
might be interesting to investigate whether second-order 
DM effects might yield an effective finite £ value. 

Using a microscopic Hamiltonian, detailed fits to the 
four magnetization step data obtained in large pulsed 
fields on powder samples of the AFM N14 tetramer 
[Mor 2 3 o(A*2-OH) 10 H 2 {Ni // (H 2 0) 3 }4], or {Ni 4 Mo 12 } 
were presented. [46J Although the molecule has C\ v sym- 
metry, it is close to exhibiting C 3v symmetry. Since the 
steps were unevenly spaced, the authors assumed the 
molecule to have weak, but important biquadratic in- 
teractions. To limit the number of fitting parameters, 
they assumed C 3v symmetry for the Heisenberg and bi- 
quadratic interactions, and Td symmetry for the single- 
ion and anisotropic exchange interactions. In addition, 
they allowed the two Heisenberg interaction strengths to 
have strong magnetic field dependencies. Subsequently, 
magnetoinfrared studies of that compound were made, 
revealing only very small differences in the responses at 
B = and 14 T,[6JJ providing little, if any justifica- 
tion for such strong magnetic field dependencies of the 
Heisenberg interaction strengths. 

More recently, a remarkably simple fit to the four level- 
crossing midpoints was made by Kostvuchenko. [50l| In 
this fit, the {N±4Moi 2 } molecule was assumed to have Td 
symmetry, so that J g = J' g = — J, and a convincing argu- 
ment was presented that the strength — J 3 of the isotropic 
three-center quartic interactions ought to be comparable 
in magnitude to that (— J 2 ) of the biquadratic interac- 
tions. Since each of these terms preserves the s quan- 
tum number, the Hamiltonian matrix in the absence of 
single-ion and anisotropic exchange interactions is block 
diagonal, and for s% = 1, it is possible to obtain the ex- 
act eigenvalues in terms of the three parameters J, J 2 , J 3 . 
However, in his fit to the magnetization level crossing 
midpoint data on {N±4Moi 2 }, he found J 2 = J 3 , which 
implies that he claimed to fit the four linear equations 
for the four level crossing midpoints with two parame- 
ters. Although two of those equations were nearly de- 
generate, three were clearly non-degenerate, rendering 
his two-parameter fit inappropriate. 

In Subsection F of the Appendix, we extended the cal- 
culation of Kostyuchenko to the five lower symmetries, 
so that there are six isotropic interactions J g , J' g , , 
for q = 1,2. In the limit investigated by 



and J 9 lq 
Kostyuchenko, J g = J' — —J, 



J 9 

b,l 



- J 9 
~ J b,2 



Kostyuchenko for the s — 4 and s — states, and one 
each of the s = 2 and s — 3 states. However, our re- 
sults do not agree with his for the other s = 1, 2, 3 states, 
except for the special case J 2 = J3, which is what he 
claimed to have obtained in his inappropriate fit.(Ho| Al- 
though Kostyuchenko presented no details of his calcula- 
tions, and provided no reference to the matrix elements, 
in Subsection F of the Appendix, we provided the explicit 
details of our calculations. We note that his fit assumed 
no widths to each level crossings. Although this is essen- 
tially correct for the first two level crossings, it is certainly 
not true for the third and fourth level crossings. 

Even if one takes the correct forms for the eigenstate 
energies with Td symmetry (neglecting the second or- 
der single-ion anisotropy contributions) given in the ap- 
pendix, one still has to solve four equations with the 
three parameters J, J 2 , and J3. In the case that there 
might be an accidental remaining degeneracy, we have 
tried to do this. We first assumed J 2 > J 3 , but no 
consistent solution to the four level-crossing equations 
could be found. We then assumed J3 > J 2 . In this 
case, the minimum energies are Eq = 16 J% — 8J3, £1-, 
E 2 = 3J+ f J 2 - f J 3 , E 3 = 6J 
E± = 10 J + 6 J 2 + 12 J 3 , where £i_ = J ■ 
±^f{§lJ 3 - 75J 2 ) 2 + 320(3J 2 - 4J 3 ) 2 . In this case, the 
square root appears in the equations for the first and sec- 
ond level crossings, so that we add those two equations 
to obtain three linear equations in the three unknowns. 
Solving these three unknowns, we then find 



X J 2 - |^3, and 



87 t 
-5-J2 



31 7 
T J 3 



J/k B = 7.511K, 
J 2 /k B = 0.7637K, 
J 3 /k B = 0.9665K. 



(79) 
(80) 
(81) 



= -J2 



and J 9 ! 



J 9 

J t,2 



-J3, our results agree with those of 



These results necessarily fit the midpoints of the third 
and fourth level crossings precisely. Now, substituting 
these values into the equations for the first and second 
level crossings, we obtain 4.35 T and 9.05 T, which are in 
remarkably good agreement with the experimental val- 
ues of 4.5 T and 8.9T, respectively. Hence, although 
Kostyuchenko didn't obtain a correct set of formulas or 
a correct fit, his idea that the three-center terms could be 
important is valid. [13] It is remarkable that one is able to 
a good fit to four level-crossings with only three param- 
eters. Nevertheless, these parameters do not give rise to 
any widths to the transitions, unlike the experiments, fill] 
In Subsection F of the Appendix, we listed the Type- 
I and Type-II first-order level crossing inductions for 
s% = 1 AFM tetramers quantization scheme is appro- 
priate for C\ v symmetry, the number of independent pa- 
rameters for that low symmetry is very large. Never- 
theless, one can quantitatively fit the four experimental 
level-crossing induction midpoints at 4.5 T, 8.9 T, 20.1 
T, and 32 T, by assuming some approximate symmetry 
such as Did or S4, for which J| is non- vanishing. The 
midpoint of the level crossings occur at 8 = tt/4, and 
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the level-crossing widths are obtained from the differ- 
ences between the values at 6 = and 9 = n/2. We 
first tried to fit the data assuming a Type-I tetramer. In 
this case, the widths of the four level crossings are de- 
termined by the single parameter | Jf — 33 9 s \, and are in 
the proportions 49 : 19 : 65 : 105, respectively. Taking 
| Jf - 3J c 9 ff | = 2.OT7, we obtain the widths to be 1.4T, 
0.54T, 1.85T, and 3.0T, which overestimates the width of 
the first level crossing, and underestimates the widths of 
the third and fourth level crossings. Nevertheless, it is in- 
teresting to try to fit the midpoint data with this assump- 
tion. We note that for Type-I, the level crossings do not 
depend upon 3 9 2 and J' g , and 3 9 2 only enters the level- 
crossing equations via the combination 3 g + 23 9 2 . We 
first assumed Jf — 3 3 9 s = +2.OT7. In this case we could 
fit the four equations with the remaining four parame- 
ters, and obtained (J g ~ 2J° 2 )/k B = -11.02K, 3 9 /k B = 
-3.522K, Jf ff /fc B = -2.069K, 3{ Jk B = -2.752K, and 
3 9 /k B = 



3^ 2 /k B 



3 9 2 /k B = 0.614 K. 



(93) 



0.503K. This fit gives a large value to the 

jotropic exchange J 
-2.OT7. This led to 

(J 9 



symmetric anisotropic exchange 3 9 s . We then assumed 

J9 — 3 J 9 
,j z oj cff 



2Jl 2 )/k B = 


-10.24K, 


(82) 


3% lk B = 


-2.570K, 


(83) 


J? S /k B = 


+0.039K, 


(84) 


3 9 Jk B = 


-1.861K, 


(85) 


JlJk B = 


-0.784K, 


(86) 



which is a much smaller value of 3 9 s = 



and a smaller ratio of 3 9 1 



to 3%,, 



which is reason- 
able. However, the fit to the four level-crossing widths is 
mediocre, at best. 

It therefore seems that with the nine parameters in the 
four Type-II tetramer level-crossing inductions listed in 
the Appendix, one might do better by fitting not only 
the midpoints but also the widths of at least three of the 
level-crossing inductions. The four widths are governed 
by the three parameters 3 9 z and Jf for q = 1,2. Setting 
Jf z — 3 2z = — Jf/2, for instance, leaves the widths in 
the proportions 3 : : 26 : 31, and their magnitudes 
are then set by | Jf | . We then fit Jf to the half- width of 
the fourth transition, which is roughly 8.0T. The fit with 
more reasonable parameter values is obtained for Jf < 0. 
To limit the remaining parameters, we arbitrarily take 
3^ 2 = 3 9 2 , and to insure Type-II behavior, choose 3' g = 
1.5J g . Then, we fit the midpoints of the four transitions 
with the remaining four parameters, and we find 



3 g /k B = 

3' g /k B = 

3 9 /k B = 

3{Jk B = 

Jii/kB 

3lJk B 



-6.799 K, 
-10.198 K, 
-4.17K, 
J£ z = 2.08 K, 
4.602 K, 
-2.529 K, 



(87) 



(89) 
(90) 
(91) 
(92) 



The resulting widths of the transitions are 0.8 T, 0, 6.7T, 
and 8.0T, respectively. These widths are in good agree- 
ment with experiment. [46| and the magnitudes of both 
Heisenberg interactions are larger than those of the other 
interactions, justifying the first-order perturbation fit. 

While this is certainly not the best fit, it is quantita- 
tively in agreement with experiment, and does not in- 
volve the assumption of strongly field-dependent Heisen- 
berg interaction strengths. [46j We emphasize that this 
fit is not optimized, as we made the arbitrary choices 
3 9 2 = J^ 2 , and 3' g = 1.5 3 g , although the only restric- 
tions on those parameters were 3 g < and 3' g ~3 g < 0. In 
addition, non- vanishing DM interactions (which vanish in 
first order) do give some additional widths to the level 
crossings, and these might provide an additional contri- 
bution to the broad third and fourth level crossings ob- 
served in experiment. The best fit to experiment may 
not be either Type-I or Type-II, but may involve a more 
complicated analysis involving other states within some 
of the constant s manifolds. 

However, with only polycrystalline data available, it 
is difficult to distinguish the different possible interac- 
tions uniquely. When single crystals of sufficient size 
for low-temperature magnetization measurements are 
made, we intend to fit the data using a more consis- 
tent set of parameters, neglecting any field dependen- 
cies, if possible. 51| To limit the number of parame- 
ters, we intend to make the assumption of C% v sym- 
metry, which will require a new quantization scheme, 
\v) = |s, m, S123, S12, si) and the appropriate single-ion 
matrix elements, which are not yet in the literature. (3ll| 

We note that our formulation of the single-ion ma- 
trix elements in terms of a pair of dimers is appli- 
cable to low-symmetry systems such as Mo^Onof^- 
OH)ioH 2 {Ni(H20) 3 }4, abbreviated as {Ni 4 Moi 2 },0 
systems such as N14 tetramers obtained from salts of 
[Ni4(H 2 0)2(PW9034)2] 10 " with C 2v symmetry, [Hi and 



4"- 



57 



the unequal-spin systems Mn^Mn"' and Ni^Mn 
58| In the first system with C\ v symmetry, one would 
expect many more single-ion, symmetric anisotropic ex- 
change, and DM interactions, making definitive fits to 
the existing powder magnetization data problematic. [46j| 
However, to improve the fits to C^ v or D3 symmetry 
systems, such as the Fe4 compound Fe4(thme)2(dpm)g 
and the Cr 1 M3 tetramer with an s = 9/2 ground 
state, 52 , 5^| would require a reformulation of the single- 
ion matrix elements as a trimer plus a monomer. [5 lj 
Even classical Heisenberg models of such systems show 
strongly different dynamics than of systems with Td 
symmetry. [U 

Unless they vanish identically, as for Td and G± v sym- 
metries, DM interactions will appear in the second-order 
eigenstate energies. This is true even when the DM in- 
teractions are site-dependent and average to zero, just 
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as for the cases of site-dependent single-ion and symmet- 
ric anisotropic exchange interactions. Hence, although 
they have been neglected in many fits to experimental 
data, they should be included in subsequent fits. They 
are most prominent for systems with lower symmetry, 
such as 54 or D 2 d, and the effects become increasingly 
strong with increasing si value. As an example, in the 
Appendix, we derived the symmetry-allowed DM inter- 
actions for the lower-symmetry C%%, appropriate for the 



[2x2] grid tetramers. [47L 14 



Finally, the DM interactions also can give rise to an 
electric polarization, and hence to multiferroic behavior. 
Our results suggest that this behavior should apply for 
tetramers with all possible individual spin values, as long 
as there is no center of inversion symmetry connecting 
the interacting spin pairs. For tetramers with S4 and 
Did or lower symmetry, this should be observable. Our 
results indicate that these effects should also occur for 
quantum spins. In addition, there is an interesting parity 
effect present in the systems we studied, with |P S (0)| cx 
(8 — 7r/2) 2si+1 . This deserves further study to elucidate 
its generality. 



IX. CONCLUSIONS 

We presented a theory of high-symmetry single 
molecule magnets, including a compact form for the ex- 
act single-spin matrix elements for four general spins. We 
used the local axial and azimuthal vector groups to con- 
struct the invariant single-ion and symmetric anisotropic 
exchange Hamiltonians, and the molecular representa- 
tion to obtain the Dzyaloshinskii-Moriya interactions, for 
equal- spin si tetramers with site point group symmetries 
Td, D^h, D2d, 5*4, Ckftj or C$ v . Each vector group intro- 
duces site-dependent molecular single-ion and anisotropic 
exchange interactions. Assuming weak effective site- 
independent single-ion, symmetric exchange anisotropy, 
and isotropic biquadratic exchange interactions, we eval- 
uated the first-order corrections to the eigenstate en- 
ergies. Depending upon the relative strengths of the 
near-neighbor and next-nearcst-ncighbor Hciscnberg ex- 
change interactions, there arc generally two types of high- 
symmetry tetramers. For the single-ion and symmetric 
exchange anisotropy interactions, we provided analytic 
results and illustrations of the antiferromagnetic level- 
crossing inductions. We also provided Hartree expres- 
sions for the magnetization, specific heat, EPR absorp- 
tion, and INS cross-section, which are accurate at low 
temperatures and arbitrary magnetic fields. For fer- 
romagnetic tetramers, we provided a procedure for a 
precise EPR determination of three of the microscopic 
anisotropy parameters. We predict that geometrically 
frustrated tetramers with symmetries S4 and Did-, as 
well as Cvj^j are likely candidate materials for multifer- 
roic states. Our procedure is extendable to more general 



systems. 
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APPENDIX 
A. Symmetry operation matrices 

Rotations by ±7r/2 about the z axis are represented by 
/ ±10 



01,2 = 



Tl 
\ 1 



(94) 



Rotations by it about the x and y axes are represented 
by 







3.4 




(95) 



Rotations by 7r about the z axis and reflections in the xy 
plane arc respectively represented by 



05,6 = 




(96) 



Rotations by 7r about the y = ±x diagonal axes are rep- 
resented by 



07,8 




(97) 



Reflections in the xz and yz mirror planes are represented 
by 



09 




(98) 



Reflections in the mirror planes containing the z axis and 
the diagonals y — ±x are represented by 



±1 
0u.i2 = I ±1 
1 



(99) 



Reflections in the mirror planes containing the y axis and 
the lines z = ±x are represented by 







13,14 




(100) 
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Reflections in the mirror planes containing the x axis and 
the lines y — ±z are represented by 



015, 



16 




(101) 



For Trf, clockwise rotations by 2w/3 about the cube diag- 
onals are represented by 



O 



17,18 




(102) 



and 



O 



19,20 




(103) 



Counterclockwise rotations by 27r/3 about the cube di- 
agonals are represented by 0\ = O^ 1 for A = 17, . . . , 20. 
Finally, there are the six 5*4 improper rotations consist- 
ing of rotations about a high- symmetry axis by ±7r/2 
followed by a reflection in the plane perpendicular to the 
rotation axis. For S 4 symmetry, z is the high symmetry 
axis, and the operations are represented by 



021, 



22 




(104) 



For Td symmetry, we also have 



O 



23,24 




(105) 



and 



O 



25,26 — 




(106) 



As described in the text, C 4 h symmetry involves 0\ = 
(0 2 ) T and 6 . 

1)4/1 symmetry contains the same three symmetry op- 
erations of Cih symmetry, rotations by ±tt/2 about the 
z axis, 01,2, and reflections in the xy plane, Oq. In ad- 
dition, it is also symmetric under the four rotations by tt 
about the x and y axes, represented by 03 a. and about 
the y = ±x diagonals, represented by 07 s.[33|| 

For C 4v symmetry, there are six group operations. 
These are rotations by ±7r/2 about the z axis, 0i,2, re- 
flections in the xz and yz planes, represented respectively 
by 09,10, and reflections in the diagonal mirror planes 



containing the z axis and the lines y = ±x, represented 
by 0ii,i2- 

For the lowest group symmetry under study, S4, the 
only two group operations are clockwise and counter- 
clockwise rotations by tt/4 about the z axis, followed by 
a reflection in the xy plane. [33I] These improper rotations 
are represented by 02i,22- 

Besides the identity operation, D 2 d group symmetry 
has five operations. The first three are rotations by 7r 
about the x, y, and z axes, respectively represented by 
03,4.6- In addition, there are two diagonal mirror planes 
associated with the principal axis, z. These are repre- 
sented by 011,12- 

Finally, the highest symmetry under study, Td, has 23 
operations besides the identity. The first five are the 
same as for D 2 d symmetry: rotations by tt about the x, 
y, and z axes, and reflections in the two diagonal mir- 
ror planes associated with the z axis. Then there are 
the four diagonal mirror planes associated with x and 
y axes, represented respectively by 0i3, 14, 15,16- Next, 
there are the four clockwise and four counterclockwise 
rotations by 27r/3 about the cube diagonals. The four 
clockwise rotations are represented by 0i7,is,i9,2o, an d 
the four counterclockwise rotations are represented re- 
spectively by 0f 7 ,i 8 ,i9,2O = °l7A9,,waO- Finall y, there 

are the six improper rotations by ±7r/2 about the x, y, 
and z axes, represented by 021,22,23,24,25,26- 



B. Molecular single-ion interactions 

The site-independent interactions in the molecular rep- 
resentation are 



jf 4 



for g = C 4 ii,D 4h , 
Jf(3cos 2 6>?-l) + 3Jfsin 2 0? 
for g = C 4v ,D 2d , 



tS 4 

-|-(3cos 2 flf 4 - 1) 



Jf 4 sin 2 6f 4 cos(2^i 



84 \ 



= 0. 



(107) 

(108) 

(109) 
(110) 



For g = D 4 h, the only non- vanishing site-dependent 
single-ion interaction is 

K%»(n) = (-1)" + Vf-. (Ill) 

For g = C 4 h, the two non- vanishing site-dependent 
single-ion interactions are 



4v 4h = J e C4fc cos(2 X f 4 "), 



(112) 



K^(n) = (-l)» +1 J e c -sin(2 X f"). (113) 
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For g — C± v , Did-, the three non- vanishing site-dependent 
single-ion interactions in Eq. (|27|) are 



K g xy (n) = 



(_!)"+! 



+Jf(l+cos 2 0?)), 



(114) 



[(2n-l)7r/4]sin(2flf)(J»- Jf), 



(115) 



tf«(n) = -sin[(2n-lW4]sin(20?)(Jf-Jf). 



(116) 



For 54, the single- ion site-dependent interactions are 



t 54 



-JiCos(2^)f 4 ) 

+J 2 sin(20f 4 ), 
K%(n) = (-l)"(jisin(2^) 

+J 2 cos(2^f 4 )), 

AS(n) = (-ir +1 [j 3 cos(0f^^) 

-J 4 sin(^f«-^ 



O) - (-1) 



n+l 



J: 



+ J4 COS 



Jl = 2( J a 4sin20 f 4 



(sin(<^ 



s 4 _ ^ 
2 7 

) 



(117) 
(118) 

(119) 

(120) 



- jf 4 (l + cos 2 flf 4 ) cos(2^f 4 )) , (121) 



J 2 = J e i4 cos6»f 4 sin(2^f 4 ), 



■h 



Isin(20f 4 )[J a S4 + Jf 4 cos(2^f 4 )], 



J 4 = Jf 4 sin 6»f 4 sin(2 tpf 4 ). 



(122) 

(123) 
(124) 



The site-dependent single-ion interactions for Td sym- 
metry are easily obtained from those Eqs. (|114p - (|116|) by 
setting 6»f 2d -> tan -1 ^) and J e Td -> 0, 



(-1)" 

3 
72 



(125) 



= -yCOs[(2n-l)7r/4]jJ<S (126) 
K^(n) = ^sin[(2n-l)7r/4]J^. (127) 



C. Molecular anisotropic exchange interactions 

We first consider the symmetric anisotropic exchange 
interactions, letting q — 1, 2 and p = g + 1. For simplicity 



of presentation, we write 



1 



J 9 — —(19 4- T 



(128) 



Then, the isotropic exchange renormalizations may be 
written as 



SJ g = 5 J' = for g = C ih , D 4h , C Av , T d 



T D 2d 
2 /iS^ 



6J Si ,SJ' Si 



(129) 
(130) 
(131) 



sin 



^[J^ + J c S4 cos(2Vf p 4 )]. (132) 



In Eq. (|132|) . <? = 1,2 corresponds to SJs 4 ,SJ' Si , respec- 
tively. 

The non-vanishing site-independent symmetric 
anisotropic exchange interactions in the molecular 
representation are 

js = -J? for g = dh,D4,h, C& v , (133) 



iD 2d 
J l.z 



jD 2d 



(1 - 3 cos 2 6»f 2 2d ) 





sin <9 12 2d , 


jD 2d _ 
2,3 — 


jD 2d 

2,- 


J e,2 ' 


754 

8i z ~ 




-3cos 2 flf p 4 






sin 2 0f 4 cos 




2 



(134) 
(135) 



For r/ = C4/1, the non- vanishing site-dependent sym- 
metric anisotropic exchange interactions in Eq. (|32p have 
strengths 

K°? y = J c c -sin(2 X f p n (137) 



4,m ~~ ^c,n h C0S (2Xlp )j 



(138) 



where x? p = 4>i p + For g = D ih ,C iv , the non- 

vanishing site-dependent symmetric anisotropic exchange 
interaction strengths are 



I 9 - -J 9 

J l,xy ~ J c,V 

K 9 - - T 9 

*^2,xy ~ J c.2- 



(139) 
(140) 



Again, the more interesting group is g = S±. We find 

Jt%y = -JlCOs(20f p 4 ) 

+J 2 sin(20f p 4 ), (141) 
K s q % y = J lS in(20f p 4 ) 

+J 2 cos(20f p 4 ), (142) 



K S q % z {n) = -J 3 cos(0f 4 -— ^ 
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- J a sin 



K I%M 
Jl 

■h 
Ji 



-J 3 sin(^-^ 
— J a cos(^ 

-J^l + cos 2 ^)cos(2^f;) 
J^cosflf'sin^f*), 



(143) 
(144) 

(145) 
(146) 

\ sin(20f;)[J / s ; i + J c % cos(2^)](147) 
sin(2^f;). (148) 



Jc,q Sm "la 



For g = D 2 d, the non- vanishing site-dependent 
anisotropic exchange interaction strengths are 



T D 2d 
J l,xy 



J 



f,l 



1 U 12 



J, 



c.l 



(1 



COS 



3-D2, 

7 12 



). 



cos(mr/2)J®l d sin26>f 2 2ti , 
- sin(n7r/2) jfi d sin 26>f 2 2 ' 



K. 



D 2d 
2,xy 



J, 



D 2d 
2,+ 



(149) 

(150) 
(151) 
(152) 



For g — Td, there are no symmetric or antisymmetric 
anisotropic exchange interactions. 

The antisymmetric anisotropic exchange interactions 
in the molecular representation are given for g — 
C 4 fe, D± h , #4, and D 2 d by 



of side a with the position vectors relative to the origin 
given by 



x(sm[(2n - l)7r/4] + cos[(2n - 1)tt/4] cos6> 



a 

+y cos[(2n - l)7r/4] sin O 



(160) 



where the acute angle #0 satisfies 7r/2 > 0q > 7r/3 for 
the [2x2] grid compounds. 4^, 5H There are three 



symmetry operations O x 2 " . The first describes rotations 



about the z axis by 7r, equivalent to 0^ 2d . The other two 
are mirror planes containing the z axis and the diagonals 
of the rhombus. [49l| which are described by 



y 2,3 



^ ± cos 6*o ± sin 9 
± sin 6*o T cos 6> 
V 1 



(161) 



For the NNN and next-next-nearest-neighbor DM inter- 
actions, corresponding to pairs across the diagonals, van- 

c 13 

ish due to Moriya rule (3) and invariance under 2 §° ■ 
However, as for dh and D^h symmetries, the DM inter- 
actions between NN spins do not vanish for C 2 % symme- 
try, but are given by 



H 



DM 



4 

E 

71=1 



\n+l 



-dcos(rt7r/2) 



z + dsin(n7r/2) 
S n x S 7l . j 



(1 = On 



(162) 
(163) 



for a general two-vector d in the xy plane. We note that 
Eq. (|162p is invariant under all three symmetries of C 2 % . 



df(n) 


= d% 


for g 




(153) 




= 


for g 




(154) 


d%{n) 




for 5 


— 5*4, -D2d, 


(155) 


df 2d 








(156) 




- + y), 






(157) 


df 4 


= d s x tx + d S y $y, 






(158) 


df 4 


= d S x 4 2 x + d s y ly. 






(159) 



Tetramers with the lowest-symmetry S'4 require five pa- 
rameters to describe the full DM interactions, those with 
D 2 d symmetry require three parameters, those with ei- 
ther dh or Dih symmetry require just one parameter, 
and tetramers with Td or Ga, v symmetry have no DM 
interactions. 



D. C\l DM interactions 

For the [2x2] grid compounds with approximate C 2 ^ 
symmetry, the four spins lie on the corners of a rhombus 



E. Compact single-ion matrix elements 

By using the Schwinger boson technique of represent- 
ing a spin by two non- interacting bosons, and check- 
ing our results using the standard Clebsch-Gordan al- 
gebra with the assistance of symbolic manipulation soft- 
ware, we find the single-spin matrix elements with gen- 
eral {s n } = (si, s 2 , s 3 , S4) to be 



(v'\S n z\v) = S„ 



m5 s ', s r 



13 ,S 13' S 24' S 24 
{s n } , — S — 1,71 



'13'°13' £, 24' £, 24 



1 A ,0 -L S C a' a a' 

a 13> a 13' 8 24' a 24 



(164) 



t{s ti } ,s,n 

S 13 ' S 13' S 24' S 24 

{s n },— s — l,rt 



13' a 13' a 24' a 24 
cr.m A {Sn},S,Tl \ 



■ 3 13' 13' 24^24 



(165) 
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C s m = Vs 2 -™ 2 , (166) 

D°> m = ay/{s-am)(s-am-l), (167) 

+^ 13 , Sl3 4<:r 4 ,n (s24,44), (168) 

$s' 24 ,s 24 e n Psust' n ( s 13i S 13) 

+^ 13 , Sl3 e+te S ' n (s 24 ,44), (169) 



E/V,S13 + (1+<X)/2,S 2 4 
V Si ,Sl,(T(T / 5+(r(l+(T / 



)/2 



a'=±l 
+ (si 3 <-> S24) 



(177) 



S 13' S 13' S 24' S 24 



a s 2 i%3 ( S 13' S 13) — T(l + £s,S 13 ,S 24 )$s' 13 ,S 13 



-V2sin[(2n - 1)tt/4] (f^S^.^ 



ct'=±1 

±8 y [fG si3+(1+ , CT) / 2 ^ , 

/ / I \ si,si,<7<7's+er(l-(-(r 
<r=±l 

+ (si 3 <-> S 24 ) ), 



2 ( 2 

7 ?S+(l + (T ')/2,Si3,S2 



)/2 



~- l Sl,S3,S %,Sl3 + l|: 

„s 24 ,s,n< / x _ (~1) r 

Psi,s 3 vSl3jSi37 — . r /s,si3 : S24 s' 13 ,s 1 3 



(170) 



(178) 



where the Ffffffi, <?£^*, Tfe,^, and £^,5,, are given 
byEqs. p72| - (fl75|) , respectively. 

In order to calculate the (v\h0l t \u) we first write it as 



+V2sin[(2n - 1)tt/4] (G^fy,,,^ 

(171) 



>\Ki\ v ) 



i 9 



(^'I'S'n ' Sn+1 + SVi+l ' S n \ 



E 1 Sl3,S 2 4 
Sl,S3,S 



/^Sl3,«24 _ 
Sl,S3,S 



Cz,x,y 



I /^fSl3 + l; s 24 X 
+ 1*31,33,-3 °3> 13 ,3 13 +1 



„ 4S24 4S-S13 

'/S13,S1,S 3 ^1- S 4-S13^ 1 S24 



(179) 



4s(s + l) 



„ 4S24 4S24 

'/si3,Si,S3^1 s + Sl3 ^l s + Sl 3 



4sa/4s 2 - 1 

^-x-\-z y 

v/z 2 (4z 2 -l) = ' 
s(x+l)-y(z/ + l) 
«(z+l) 



;[1±(-1)"], 



(172) 

(173) 
(174) 

(175) 
(176) 



where A™ is given by Eq. (|46|) . The prefactors m, AJ m , 
C s m , C™.^, Df m , and D*'™! are consequences of the 
Wigner-Eckart theorem for a vector operator. [33} The 

challenge was to obtain the coefficients r^ s "V'™ , and 

S 13' S 13' S 24' S 24 

A^ S "V'" / . Their hierarchical structure based upon 

S 13' S 13' S 24' 5 24 

the unequal-spin dimer suggests that analogous coeffi- 
cients with n > 4 may be obtainable. [f| Details will be 
presented elsewhere. [5 lj 



F. First-order eigenstate energy constants 

The constants appearing in the first-order eigenstate 
energies (|5"Tj) are 



= - ( 1 + f 2 -r? -ii 2 

-V ^S,S13,S24 'S,S13,S24 '/S+1,S13 ,S2 



We note that for general n, S n ■ S n ±± commutes with 
the isotropic part of the Hamiltonian, which includes the 
Heisenberg, biquadratic, and isotropic three-center inter- 
actions. Thus, we expect (v'\S n • S„+i + S n +i ■ S n \u) 
to vanish unless m! = m and s' = s. It is easy 
to see that ml = m by inspection. The matrix ele- 
ments for s' = s ± 1 and s' = s ± 2 can then eas- 
ily be shown from their s,m dependencies to vanish. 
For example, the s' = s — 2 terms are proportional to 
C™C™_ ! + \ J2a D^ m D^~ a = 0. There are three 
non-vanishing s' = s terms, corresponding to the inter- 
mediate states s" — s, s + 1, and s — 1. These are respec- 
tively proportional to m 2 + \ J2 a Aj* m A i ^ m+ * ) = s(s + 
1), C™ s -iC? +1 = («+l)(2 S + 3), 

and CTC m s - \ £~ D-*> m D*:™- 5 = s(2s - 1). Hence, 
the Wigner-Eckart theorem guarantees that these are in- 
dependent of m. We first performed two checks of our 
matrix element forms. First, we evaluated 



(i/'ISi-Ssk) 



s(s + 1) 



j 2 



_( p3x3,S2i\ _ ( jps\3+l,32i\ 
1 1 S\,S\,S J \ 1 S\,S\,S j 

+ S ( 2s - 1 )[jQVl sl3 , S24 

_ I /~IS 13 ,S24 \ 2 _ f (~H> 13 + 1,S2 4 \ 

\ si,si,—s J \y r s 1 ,s 1 ,—s J 

+ ( S + l)(2 s + 3)[lr ? 2 +Mi3iS; 

_ ( S-ISl3,S24 \ 2 _ ( (-1313+1,324, 

^31,31,3+1) a u »u»+l 
^Sl3(si3 + 1) - S\{S\ + 1)^5^,^', 



( 



0± = CI. 



±±2(J2[(f:;:S +1)/2 > s 



<T=±1 



5 v y 



(180) 
(181) 



24 



as required. Similarly, {i > '\S2 ■ S^\v) is found from the 
above by setting si 3 <-> s 2 4, as required. Then, we found 

Y,{ASn-S n+1 \v) = ^5 v y, (182) 



n=l 



where 



+ l)(2 S + 3)^ +1 

>S13>S 



= 2 



s(s + 1) - S1 3 (S1 3 + 1) - S24(«24 + 1) 



(183) 



as required. We then may write 



(u'\S n ■ S n+1 + S n+1 ■ S n \v) = 5 B , a ,5 m , m >M™g(V, n), 

(184) 



where 



M S24 ' s , 24 (77,n) = 5„> s 5 a > . 

«13)«i3 *13' s 13 s 24>°24 



9o(v) 



+ E 5s 24. S 24+ ff ^i3. S 13+ ff ' 

<7,cr'=±l 

-V2sin[(2n - 1)tt/4] 



V <5 S ' s 5.' s 

L / s 13' s 13 s 24' s 



9<y{v) 



CT=±1 



-V / 2sin[(2n + 1)tt/4] 



E* 



T<5.«' 



CT=±1 



13' S 13+ <T °24'°24 4 



(185) 



where 



S(S + 1)(1 + £ S;Sl3iS24+CT (i +(7l )/2) 



E 

<7l=±l 



x pS24 + (l+c)/2,Sl3 
Sl,Sl,S 



+ i y (2s + 1 + cti)(2s + 1 + 2(7!) 

CT2=±1 

X % + (cri+l)/2,Sl3,S24 + (CT + <T 2 )/2 



xG 



S 2 4 + (l + ff)/2,Si3 
Sl,Si,cri<7 2 s+cr 2 (l+CTi)/2' 



(186) 



CT = s(s + l)^ Sl3 + ( CT + 1 )/ 2 ' S24 i^ S24 + ( 1 +' T )/ 2 ' Sl3 + <T 

+ - J2 (2s + l + cn)(2s + l + 2(Ti) 

<Tl=±l 

X (7 ,Sl3 + ( 1 +' T )/ 2 ' S24 /^rS24 + (l+CT)/2,Sl3+(T 

Sl,Sl,0 , i(TS-fO'(l + (Ti)/2 Sl,5l,— C7l(TS — (j(l+(Tl )/2 

+ (*13 ^ *24), (187) 

/*+,_ = £ s( S + l)^ s -+i' SM - (1 - CT)/2 F^+( 1 +-)/ 2 

CT=±1 



+ i (2s + cti + 1)(2s + 2c7i + 1) 

CTi=±l 

x g,S13 + l,S24-(l-<T)/2 £jS 2 4,S13 + (l + <x)/2 

Si ,Sl,CT(TlS+(7(l+<Tl )/2 Si ,Si ,(7(TlS+(T(l+Cri )/2 ' 

(188) 

and where g<j(p) and /i_ j+ (z7) = h+-(p) are respectively 
obtained from .g<x(z7) and the /i +j _(z7) by setting si 3 <-» 

S24- 

Letting a; = s i3 and y = ,s 24 , these expressions may be 
simplified to yield 



h a ^(0,X,x) = -y/(2x + 2 + d){2x + (T)vl+( 1+<T )/2,suSl 

^ 2x + 1 + °l, (189) 



2 

(x - vf 

s(s + l) 



Vx+(l+a)/2,s 1 ,s 1 Vy+(l+a)/2,s 1 ,s 1 



h+,-(s,x,y) 



X\/{x + y + l + (T) 2 -s 2 
x^J(x + y + l + a) 2 -(.s + 1) 2 , 

fors>l, (190) 
{x + y + lf 



s{s+l) 
x^is + 1) 2 -(x-y-1) 
x^/s 2 -(x-y-1) 

g a {s,X,y) = Vx+(l+a)/2,s uSl 



?7x,Sl,Sl?7l/+l,Sl,Sl 

2 

(191) 



x^[x+(l + a)/2] 2 ~{y-sf 
x v /(y + .s + l) 2 -[ 2 ;+(l + ( T)/2] 2 . 

(192) 

The diagonal matrix elements of l are then easily 
found to be 



(193) 



where 



(7=±1 



+h\JV) + h 2 +> _{V) + -g 2 (V). (194) 

The three-center isotropic quartic spin-spin interac- 
tions proposed by Kostyuchenko may be written for the 
six g symmetries as 



(195) 



Wf l — — g^l E {Sn- S n+ \+ Sn+1- S r }j 

n=l,3 

X ^] ^S„/ • 5„/ + i + S n i+\ ■ S n '^j 

n'=2,4 

+ ff.C., (196) 



W? 2 



— — - J t 9 2 ^1 • S 3 + S 2 • & 
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5 


16 

5 
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45 


43 
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1 
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11 
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3 
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IV: Values of 


c_ 

V 


, #77, 


and 7j7 for si = 1/2. 
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2 
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1 
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3 


(i 


6 
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4 
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1 


2 


2 


7 
5 


14 

5 


36 
5 


32 
5 


21 
10 


21 
10 


31 

4 


2.', 
4 






x] 

n 


=1 






+1 ' 


+■ s„ 


+ 1 ' 
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1 
1 


2 
1 


1 
1 


3 
5 

1 


8 
5 

-2 


(18 
15 

4 


32 
15 




11 
10 
1 
2 


9 
10 
1 
2 


23 

4 

15 

4 


1 

12 
1 

4 



where for g = Td, we have jj\ = = — J3 in the 
notation of Kostyuchenko.JE^ From our matrix elements 
above, it is then easy to see that 

/ >\u9 1 \ 9o\p)Sv,v' T g ( , . ,x 

\ V \m,%\ v ) = ^-^^ t2 ( s 13(si3 + l) 

+*m(«m + 1)-4si(si + 1)), (198) 



1 1 
2 
1 







_ 11 

3 
U 
3 



2 

6 

22 
3 

_ 8 
3 



2(1 

3 

16 
3 

16 
3 

16 
3 



-20 

3 

8 
3 

_ 8 
3 

_ 16 

3 



1 

2 

9 
2 

11 

(i 
1 
2 



9 -# 



2 

11 
6 
1 
2 



32 
3 



ill 

3 



23 
3 
22 
3 



16 
3 



TABLE V: Values of cd 



*)-, c-, Bj7 and T- for Si = 1. 



where go(v) is given by Eq. (|183|) . With regard to the 
matrix elements (z/j'Tij'-jJ^), we first note that 



sin[(2n ± 1)tt/4] = ^ sin[(2n ± 1)tt/4] = 0, 

(199) 



n=2,4 



n=l,3 



so that there are no contributions from g a and g a . The 
diagonal matrix elements (i>\Ti. 9 ^v) is then easily found 
to be 



- J 9 T- 



64 



cr=±l 



(200) 



(201) 



For si = 1/2, the matrix (u'\H^ i + Wf x \u) is diagonal, 
so these interactions can be treated exactly. For S\ > 1, 
since these quartic interactions preserve s, the resul ting 
matrix is block diagonal, as noted by Kostvuchenko. [50j 
To the extent that the single-ion and exchange anisotropy 
interactions can be neglected or treated in first order only, 
the matrix of the resulting Hamiltonian is block diago- 
nal. Kostyuchenko compiled a table of the diagonalized 
eigenstates for Td symmetry with s\ = 1, neglecting the 
anisotropy interactions. In Table IV, we calculated the 
exact eigenstates for s\ — 1 of the above Hamiltonian for 
the six symmetries under consideration. 

In Tables IV- VII, and VII-XI, we have listed analytic 



formulas for the coefficients c£, b£, and ci, and B- for 



the lowest four eigenstate manifolds of FM and AFM 
tetramers, respectively. In both cases, the manifolds are 



V = SS13S24 




422 


-4J 9 - 6J' g - 4J 9 , - 2J 9 2 - 4J 9 X - 8J 9 2 


322 


-6^-4J b 9 iX -2J b 9 i2 


312,321 


9 7 AT' 12. T 9 op ±1/9 

£Jg g J h l ^b,2 ' 9 J t,l 


212,221 


f aJi 307 t9 nrS 1 55 79 
,J 9 36 J h,l b,2 36 t,l 


211 


7 9 7' 31 79 n 79 i 23 79 1 9 79 

Jg £J g 4 J b l ^ b ,2 4 J t,l LJ t.-. 


222,202,220 


-3^-4^,-5^%, E 2± 


122 


2 7 — 3 f — — I 9 — 2 J 9 —-i/9 1 1(179 
LJ g OJg 4 J h l ^ b ,2 4 J t,l ' ±UJ t, 


111 


7 ni 15 79 9 79 1 79 9 79 
,J 9 £ ">g 4 J b,l i,y b,2 4 t,l ^ u t,2 


121,101 


E 1± 


112,110 


Ei± 


011 


-2J 9 + 2J' g - 8J b 9 ! - 2J b 9 2 + - 4J* a 


022,000 


£o± 



TABLE VI: Eigenstate energies £xr for si = 1 as a func- 
tion of the quantum numbers u — s,si3,S24 in the ab- 
sence of anisotropy interactions. The two additional s — 2 
eigenstate energies E 2 ± = \[a ± «/6 2 + 56(J bl ) 2 ], where 



o — ^ 7" Q 7 ; 63 79 _ 7 79 I 15 79 _i_ C 7 

u — OJ 9 »J 9 4 J bjl ' ^b,2 ' 4 u t,l < ""'l 



t 9 2 and b = 



3J 9 3J 9 + 4 J b 
of doubly-degenerate s 



+ 3JZ a 



^ Jt 9 , + 6J t 9 2 . The two sets 



1 eigenstates have E±± = i[A + 
B ± v/(A- fi) 2 +4C 2 ], where A = 3J 9 



9 79 I J_ 79 

^ J b,2 ~ 12 J t,U 

^r( 3J b,i - A 79 

where a = 



q /9 



AT' _ 23 79 
^ J 9 4 J b,l 

5J 9 2 + f J 9 !, and C 



16 f 79 
3 *. J t,l 



4J 9 !), and £ ± = U a 
-6J S + 6J 9 - f J 9 4 - 
8 J 9 



+ /3± V(a-/3)2 +4 7 2], 



2<2 



f ^i + 12J t 9 2 



7 9 

b,l 



and 7 : 



4-v/S / 79 _ 79 i 
3 Wt,l J b,l) 
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S, S13, S24 



S, S13, S24 



6± 



6^ 



4si,2si,2si 
4si — 1, 2si, 2si 
4s i - l,2si,2si - 1 
4si - 2,2si,2si 

4s i - 2,2si,2si - 1 



4s i - 2,2si - l,2si - 1 
4si - 3,2si,2si 
4si - 3,2si,2si - 1 



4si - 3,2si,2si - 2 

4si - 3,2si,2si - 3 
4si - 3,2si - l,2si - 1 
4si - 3,2si - l,2si - 2 



2 S1 -1 
8s 1 -1 

(2 31 -l)(4si-3) 
(4 sl -l)(8si-3) 

(2 31 -l)(4 31 -3) 
(4 31 -l)(8si-3) 
(2si-l)ai(si) 
(4si-l)(8si-l) 

X (8*1-5) 

"2(si) 
(2si-l)(4si-l) 

X (8si-5) 

"2<£l) 
(2 Sl -l)(4 Sl -l) 

(8si-5) 

2ai-3 
8si-5 

(2si-l)(8si-15) 
(8si-3)(8si-7) 

_ "4 (si) 
(4 Sl -3P(2 Sl -l) 

* (4si-l)(8si-3) 

X (Sn'-T) 

°e(si) 
(2si-l)(4 sl -3) 



2si 
8si-l 
2gl (4sj -3) 
(4si-l)(8si-3) 
2si 
8 S1 -3 

2siai(si) 
(4si-l)(8 S1 -l) 



A (8^1-5) 
4s 1 (si-l) 
(2si-l)(8si-5) 



2^3 (si) 
(2 Sl -l)(4 Sl -l) 

1 



(8si-5) 
2si 



X ■ 



' (8si-7) 
(23i-l)a 8 ( 3 l) 
(4si-3) a (8si-7) 

(2si-3)((4 31 -5) 
(4si-3)(8si-7) 

(4si-l)(4si-3) 
i 

X 



(8si-7) 



8si-5 
2si (8si-15) 
(8si-3)(8si-7) 

4£l05(£l) 

(4 Sl -3) a (2 Sl -l) 

x 1 

(4si-l)(8si-3) 

X (8 Si 1 -7) 

a 7(«l) 
(2si-l)(4si-3) 

1 

8si-7) 

2a9( s l) 
(4si-3)^(8si-7) 

2si(4si-5) 
(4si-3)(8si-7) 

2nil(si) 
(4si-l)(4si-3) 

X (8si-7) 



x- 



TABLE VII: Values of a— for the ground and first three ex- 
cited state manifolds for FM tetramers (or the highest four 
excited state manifolds of AFM tetramers), where ai(x) = 
32a; 2 - 44x + 3, a 2 (x) = 16x 3 - 36a; 2 + 26x - 3, and a 3 (x) = 



8X* 



6x 2 +2x-l, o,a{x) = 2048x b - 9472a; J + 17344x 4 



15440a; 3 + 6924a; 2 - 1530a; + 135, a 5 (x) = 512a; 5 - 1856a; 4 + 
2416a; 3 - 1420a: 2 + 399a; - 45 a 6 (x) = 16a: 3 - 52a; 2 + 58x - 9, 
a 7 (x) = 4x 3 - 7x 2 + 3x - 1, o 8 (a;) = 16a,- 2 - 48a: + 45, 
ag(x) = 16x 3 -24a- 2 + 15x-9, a w (x) = 32a; 3 - 96a- 2 + 70a; - 9, 
and an (a-) = 16a- 3 — 16a- 2 + 5a; — 2. 



restricted by < si3,S24 < 2si and |si3 — S24I < s < 
S13 + S24- In addition, the coefficients are symmetric un- 
der S13 <-> S24. Hence, for si = 1/2, the five distinct 
allowed (s,si3,S24) states are (0,0,0), (0,1,1), (1,1,0), 

(1.1.1) , and (2,1,1). For s\ = 1, the fourteen distinct 
allowed (s,Si3,S24) states are (0,0,0), (0,1,1), (0,2,2), 
(1,1,0), (1,1,1), (1,2,1), (1,2,2), (2,1,1), (2,2,0), (2,2,1), 

(2.2.2) , (3,2,1), (3,2,2), and (4,2,2). From these tables, 
the coefficients , 6^, and ci, and £?- for all of the al- 
lowed eigenstates of tetramers with si < 3/2 are given. 
For si = 2, the values for the nine states with s = 4 can- 
not be obtained from these formulas, but the values for 
the other 46 eigenstates with < s < 3 and 4 < s < 8 
are given. 



4si,2si,2si 
4si - 1, 2si,2si 
4si - l,2si,2si - 1 

4si -2,2si,2s! 
4si -2,2si,2si - 1 
4si -2,2si,2si -2 
4si - 2, 2si - 1, 2si - 1 
4si -3,2si,2si 
4si -3,2si,2si - 1 

4si -3,2si,2si -2 
4si -3,2si,2si -3 
4si -3,2si - l,2si - 1 
4si - 3, 2si - 1, 2si - 2 



24s( 
8si-l 

8si(5si-2) 

8si-3 
8si (5si -2) 

8si-3 

4(112sJ -102s J +22si -1) 
(8si-l)(8si-5) 
4(14s^-12si+l) 

8si-5 
4(14sJ-12si+l) 

8si-5 
4(14sj-10si+l) 

8si-5 

12(48sf -74s^+34si -5) 

(8si-3)(8si-7) 

4i»i(*l) _ 
(4si-l)(4 sl -3) 

1 



8ai 
8s? 



8si-3 
8si(3si-l) 
8si-3 
8si(24sj-15si+l) 
(8si-l)(8si-5) 



8s 



(8si-3)(8si-7) 
12(6s^-8si+l) 
8si-7 
12(24sj -50sj+30si-5) 
(4si-3)(8si-7) 
4(18si62-22si+5) 
8si-7 
4(72s^-98s5+34si-3) 
(4si-l)(8si-7) 



_ l("l-l) 
8si-5 

4(10sj-8si+l) 
8si-5 

8si(5si-3) 
8si-5 

8si(40s^-51si+15) 
(8si-3)(8si-7) 

8sib 2 (si) 
(4si-l)(4si-3) 

x 1 

(8si-3)(8si-7) 

4(6s^-8si+l) 
8si-7 

46 3 (3i) 
(4si-3)(8si-7) 

24si(si-l) 
8si-7 

4b4j£l) 
(4si-l)(8si-7) 



TABLE VIII: Values of b- for the ground and first three 
excited state manifolds for FM tetramers (or the high- 
est four excited state manifolds of AFM tetramers), where 



6i(a;) = 2304a; b 
& 2 (a;) = 1024a; 4 
56a; 3 



6112a; 4 + 5968a; 3 - 
2304a- 3 + 1808x 2 



2650a^ + 552a; - 45, 
576a; + 63, 63 (a;) = 



114a: 2 + 72a; - 15, b 4 (x) = 56a- 3 - 70a; 2 + 20a; - 1. 



For Types I and II tetramers, the sth AFM level- 
crossing induction in the first-order approximation may 
be written as 

7^ (1) (9) = -J g s-Q(J' g -J g )2s S \Jl 2 
+2 ai ( fll + l)J» a )£?(!±i) 



- 1 9 

J h.2 



E 



J? 



2 
I 9 

J l,z 



-(a+ + 2b + + a+ cos 2 9) 



[c 2 - + -(b + + b-) + c^ cos 2 9)] 



js 

-—^-{02 + b~ + a{ cos 2 9) 



(202) 



O(s) is the standard Heaviside step function, E(x) is the 
largest integer in x and the level-crossing parameters , 
6 ± , and cj for j = 1, 2 and d are functions of s, si and the 
tetramer type. For Type II, the functions are different 
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S, S13, S24 



S, S13, S24 



4si, 2si, 2si 
4si - 1, 2si,2si 

4si - l,2si,2si - 1 
4si -2,2si,2si 

4si -2,2si,2si - 1 

4si -2,2si,2si -2 
4si -2,2si - l,2si - 1 
4si - 3, 2si,2si 
4si -3,2si,2si - 1 

4si -3,2si,2si -2 

4si -3,2si,2si -3 
4si -3,2si - l,2si - 1 
4si -3, 2si - l,2si - 2 



4ai-l 
2(8si-l) 

4si-3 
2(8si-3) 

16sj-12s!+3 
2(4si-l)(8si-3) 

32sj-44si+3 
2(8s 1 -l)(8s 1 -5) 

ci (si) 
2(2si-l)(4si-l) 

x 1 
X (8^1-5) 

8s 2 -10si+5 
2(2si-l)(8si-5) 

43i-3 
2(8si-5) 

(43i-l)(83i-15) 
2(8si-3)(8si-7) 

C3( 3 l) 
2(2si-l)(4si-3) 

x 1 

(83i-3)(8si-7) 

(431-1) 
2(2si-l)(4si-3) 

(83^-1831+3) 
X (S3i-7) 

16s 2 -28si+21 
2(43i-3)(83i-7) 



4sj -1 
2(83i-l) 

4si-3 
2(83i-3) 

831 (23 i-l) 2 
(4si-l)^(83i-3) 

32s 2 -443i+3 
2(8si-l)(8si-5) 

C2 (s 1 ) 

2(231-1)^(431-1)^ 



X ■ 



1 



' (8si-5) 
23i(43 2 -73i+3) 
(23i-l)^(83i-5) 

43i-3 
2(83i-5) 

(4si-l)(8si-15) 
2(8si-3)(8si-7) 

C4(»l) 

2(23i-l)^(4si-3) 2 

v 1 

(83i-3)(83i-7) 

2(3i-l)(43i-l) 
(2si-lp(4si-3) 2 
(1633-403 2 +293i-8) 



163l(3l-l)(23l-3) 

(4si-3) 2 (83i-7) 



163i— 28si+13 
2(43i-3)(83i-7) 



2(83i-7) 

16(3i-l) 2 (23i-l) 
(4si-3) 2 (83i-7) 



TABLE IX: Values of ci for the ground and first three excited 
state manifolds for FM tetramers (or the highest four excited 
state manifolds of AFM tetramers), where ci(x) = 32a: 3 — 
56a- 2 + 30a; - 3 and c 2 (x) = 256a- 5 - 640a; 4 + 576a; 3 - 240a; 2 + 
48a; - 3, c 3 (x) = 256a; 4 - 800a; 3 + 840a; 2 - 330a- + 45, c 4 (x) = 
2048a; 6 - 8960a; 5 + 15232a; 4 - 13120x 3 + 6096a; 2 - 1440a; + 135. 



for even and odd s. 

The AFM level-crossing parameters are denned accord- 
ing to 



*(2s - l)< s 'it, 



± - „«i.± 



a 9 = sa 



-( S -l)(2 S -3) a ^ s , , 

,81, ± 



1,81, ± 
S,S 13 ,S 2 4. 



- (s - l)o 8l, f , 

S 1 ,S -. o , S . 



S (2 S -l)c^- S24 
- (s -l)(2 S -3)C 1 - u , a , , 



(*-i)< 



e 



T- — T-> 



where 



{• s 1j s 13i s 24j s i} 



(203) 
(204) 
(205) 

(206) 
(207) 
(208) 
(209) 

(210) 



and the s' 13 ,s 2 4 values depend upon the tetramer type. 
In the next two sections, we evaluate the af, b^ 1 and 



Q x x "1(3;) [03 (x) -4si (£1 + 1]] 



1, a;, x — 1 



3(2i+3)(2i-l) 

-a 2 (x) + 43!(3l+l) 



5i(3;)[a; + a: 2 +43i(si+l)] 
3(2x+3)(2x-l) 

3 + 3 2 +43i(si+l) 



a 3 (3:)-43i(3i + l)(l + 23: 2 ) 4x 2 +2x 4 +4s 1 (s 1 +l)(l + 2x 2 ) 
~ ~ 5(43 2 -l) 

a 4 (j;)[a; 2 + 3:+43i(si + l)] 
21(23+3)(23-l) 

a 7 (3;)+43i(3i + l)(5-2a; 2 ) 
21(43 2 -1) 

4-3 + 3 2 +4si(si+l) 
21 

a 8 (3)[3 + 3 2 +4si(s 1 + l)] 
45(23+3)(23-l) 

an(j:) + 83i(3i+l)aio(3:) 
30(43^-l)(43^-9) 



9 5 4 (3;)[a 5 (3:)+43i(si+l)] 
21(23+3X23-1) 

o _ _ _ I a 6 ( 3: )-43i(3i+l)(5-2 3; 2 ) 

1 21(4x 2 -l) 

2 x a; - 2 3(-3+3+3 2 )+4si(3i+i) 

q „ „ S 8 (a:) [-3a 2 (:c) + 43 1(31 + 1)] 

°' 45(23 + 31(23-11 



o _ 1 a 9 (3:) + 83i(3i+l)aio(3:) 

- 1 - 30(43^-1X43^-9) 



3, a;, a; 2 
3, a;, a; 3 



^12(3)— 831(31+1)3(3: — 2) 
18(23-1X23-3) 



ai 3 (3) + 83l(3l+l)3(3-2) 

18(23-l)(23-3) 



TABLE X: Values of ajr for the ground and first three ex- 
cited state manifolds for AFM tetramers (or the highest four 
excited state manifolds of FM tetramers), where x repre- 
sents any value of S13 that satisfies < si3,S24 < 2si and 
S13 — S24 < s < S13 + S24, and 5i(x) = 3 + 4a; + 4a; 2 , 



cb(a:) = 3(— 1 + x - 



(2x + 5) (2a; - 3), 05 (a;) = 3(1 - 



&3(x) = — 3 + 6a- + 6a; , 04(3;) 



2 ), a 6 (x) = 3(5 - 



16a- 2 + 2a; 4 ), a 7 {x) = 2a; 2 (7 - a; 2 ), a 8 (x) = -33 + 4x + 4a- 2 , 
a 9 (a;) = 198 - 729a- 2 + 330a- 4 -24a- 6 , a w (x) = 33 - 37a; 2 + 4a- 4 , 
a n (a-) = 169a- 2 - 102a; 4 , a 12 (x) = 27-78a; + 63a; 2 -24a- 3 + 6x 4 , 
and 013 (a;) = 9 - 30a- + 23a; 2 - 8a- 3 + 2a; 4 . 



cf for Type I and Type II AFM tetramers. The NN 
biquadratic exchange level-crossing parameter d has dif- 
ferently complicated forms for Type I and Type II AFM 
tetramers, the general forms of which arc not given for 
brevity. 

For si = 1/2, the Types I and II first-order level- 
crossing inductions jBy^\6) are given in the text. In 
that simple example, there are no effects of single-ion 
anisotropy. Hence to illustrate the full dependencies on 
all of the microscopic parameters, we list the s\ = 1 first- 
order level-crossing inductions. For Type I, we have 



7 < c «(0) = -J g 



-JS 



(I 



1 

12 



(35J b s 1 + 13J t 9 4 ) - 2J t 9 2 



10 



- cos 2 6 



7< C(1) W 



+ -J c s ff (l-3cos 2 tf), (211) 
-2J 9 + 2^ 1 + Jf(|-^cos 2 0) 
+4^ 1 -4J« 2 -^^ ff (l-3cos 2 0), 



(212) 



7 < c(1) (tf) = -3J g + 1 (7Jl 1 -9Jl 1 )-6Jt 2 
TO /181 13 2n 



28 



S, Sl3, S24 



b± 



b- 



0, x, x 




8si(si+l) 
3 


4i(a; + l)-8si(si + l) 
3 


l,x,x 




4[-l + a;+a: 2 +2s 1 (si + l)] 


8\x(x+l)- 31 ( 31 + l)] 




5 


5 


l,x,x- 


i 


4[b 1 (:r)-2s 1 (s 1 +l)(l-8 :I: 2 )] 


4[6 2 (a:) + 2s 1 (s 1 +l)(l-8a: 2 )] 


5(4x^-1) 


5(4x 2 -l) 


2,x,x 




b 3 (x)+8si(si+l)6 4 (x) 


b 5 (x)-8s 1 (s 1 +l)b 4 (x) 




21(2a: + 3)(2a:-l) 


21(2i + 3)(2a:-l) 


2,x,x- 


i 


4[S 6 ( a: ) + 2 sl ( sl +l)(l+8 a: 2 )] 


4[5 7 (a;) + 2s 1 ( S1 +l)(l + 8:r 2 )] 




7(4x^-1) 


7(4x^-1) 


2,x,x- 


2 


4[-3+a:-a; 2 +6si(si + l)] 


4[1-2i+2i 2 -6si(si+1)] 


7 


7 


^ ^ X ^ X 




8lbg(x) + s 1 (s 1 +l)b g (x)] 


4[6i O (a:)+2si(«i+l)Ma0] 




45(2x+3)(2x-l) 


4B(2x+3)(2x-l) 


3, x 5 x 


1 


4[6 11 (a ; )-2 si ( 31 +l)6 12 (x)] 


8lb 13 (x)+s 1 ( S1 +l)b 12 (x)] 


15(4a; :! -l)(4i ;4 -9) 




3, x ; x 


2 


4[-3+2si(si+l)] 
3 


4[a:(l-a:)+2si(si+l)] 
3 


3, x , x 


3 


4[6i 4 (a:)+2si(si+l)6iB(!c)] 


8[b 16 (x)+ S1 ( S1 +l)b 15 (_x)] 


9(2x-3)(2x-l) 


9(2x-3)(2x-l) 



TABLE XI: Values of b— for the ground and first three 
excited state manifolds for AFM tetramers (or the high- 
est four excited state manifolds of FM tetramers), where 
bi(x) = l-2x 2 - 2x 4 , b 2 (x) = -3x 2 (l - 2a; 2 ), 6 3 (a;)12(15 - 
19a; - 15a; 2 + 8a; 3 + 4a; 4 ), b 4 (x) = 9 + 20a; + 20a; 2 , b 5 (x) = 
-16x(9 + x- 16a- 2 - 8a; 3 ), b 6 (x) = 5- 16a- 2 + 2a; 4 , b 7 (x) = 
7x 2 - 10a; 4 , b 8 {x) = 3(33 - 37a; - 33a; 2 + 8a- 3 + 4a- 4 ), b a (x) = 
87 + 44a,- + 44a; 2 , b w (x) = x(lll + 43a; - 136a; 2 - 68a; 3 ), 
bu(x) = 3(-66 + 243a; 2 - 110a; 4 + 8a; 6 ), b 12 (x) = 87 + 
52a,- 2 - 64a; 4 , b 13 (x) = a; 2 (107 - 151a,- 2 + 44a- 4 ), b 14 (x) = 
-3(9 - 26a; + 21a- 2 - 8a; 3 + 2a- 4 ), b w (x) = 9 - 32a; + 16a; 2 , 
and b 16 (x) = -9 + 30a,- - 35a,- 2 + 20a; 3 - 5a- 4 . 



S,S13, 


S'2 


4 






0,x, 


X 




3+4x+4x A 
6 


3+4a:+4a: 2 
6 


l,x, 


X 




3-4a;-4a: 2 


3-4:c — 4a; 2 




6 


6 


1, X, X 




1 


3+2a; 2 


3(1-* 2 ) 




10 


10 


2,x, 






15 — 4a:— 4x 2 


15 — 4x— 4x 2 


X 




42 


42 






1 


15~2x 2 
42 


15 + a: 2 
42 






2 


13-4a:+4a: 2 


8(2+x-x 2 ) 




42 


63 


3, x, 


X 




33-4a;-4i 2 


33-4i-4a; 2 




90 


90 


3 , x , x 




1 


ll-z 2 
30 


132-17X: 2 
360 


3 , x , x 




2 


1 

3 


23+4a:-4a; 2 
72 


3 , x , x 




3 


ti-2x + x 2 


5(3+2a;-:c 2 ) 




18 


72 



TABLE XII: Values of c£ for the ground and first three ex- 
cited state manifolds for AFM tetramers (or the highest four 
excited state manifolds of FM tetramers) . 



-§^ ff (l-3cos 2 0), 



(213) 



7 < c(1) (0) = -4J g -4^ 1 + J|(|-|cos^) 

-8^ 2 -^e ff ff (l-3cos 2 0). (214) 
The Type II first-order level-crossing inductions for si = 



1 are 

7 < c(1) (0) 



-Jg- o(H<l + 9J t 9 l)+3J b S 2 



-^f (5 - 3 cos 2 0) + -ffl(l + 3 cos 2 6>), 
6 6 



(215) 



7 < c(1) (0) = -^-4 + ^(15^ + 23^ 



12 

+ ^(5 + cos 2 0) + '^±(1 + cos 2 6) 

J9 

+3J£ 2 + 2J t 9 2 + -|^(1 + 3 cos 2 0), (216) 



jB^\e) = -^-2^ + ^(107^ + 179^) 

+^(|-^cos 2 0)+^(l + 5cos 2 0) 

(217) 



7 < c(1) 



(<?) 



-2,7^ + ^ (l+cos 2 ^ 

2J 3 - 2J' g + 1(77^-29^,0 - 8J t 9 2 
83 3 9 ^\ 2 



+J '(|-r os20 ) + ^ (1+5cos20) 

(218) 



J9 

+-|^(1 + cos 2 9) 



G. Type I First-order AFM level-crossing constants 

The Type I constants are relevant for both J' g — J g > 
AFM level-crossing inductions and for some low energy 
FM manifold states (with J' g > J g > 0). The low-energy 
states within an arbitrary s manifold have S13 = S24 = 
2s\. We let fflg*^ S24 = a^, etc. For general s, si, we 
have 



,«i,± 

i s,2s i ,2si 
J s,2s 1 ,2s 1 



- <2 Sl , 2si (l + ^7— y) 



(219) 



1 



s(s + 1) 



"g,2si,2si 



K?Si 

s,2si,2si 



2(2s + 3)(2s- 1) 
+[2s(s + 1) - l][8«i(2si + 1) - a(a + 1)] 

± 4^T ( 2[16fi2 + S(S + 1)1 [S{S + 1} " 1] 
-8«i[2s(s + l)-l]) 
3[s(a + l)-l]-8ai(2«i + l) 



16 



2(2s- l)(2s + 3) 
[s(s + l)-4s 1 (2s 1 + l)] 2 



(220) 
(221) 



+ 



+ 



s(s + l)(4s 1 -s)(s + 4s 1 + l) 



8(4«i - 1) 



s 2 (4si + l) 
(4^1 - 1) 



(222) 



29 



s\{A Sl + 1) 
' (4si-l) 



$s,0, 



(223) 



From these expressions, we may evaluate the Type I 
first-order level-crossing inductions for AFM tetramers. 
From the definitions of the level-crossing constants in 
Eqs. (|203[) - (|207[) . we rewrite them to explicitly indicate 
the Si, s and type dependencies, and for Type I, it is easy 
to show that 

a?A s ) = ^"W( 1= f=4^)" (224) 
2s[8si(2 Sl + 1) + 4s 4 - 10s 2 + 3] 



x it 



(4s 2 -l)(4s 2 -9) 
1 



C ?,2 ( S ) 



4s i - 1 

3[4s 3 + 5s 2 - 3s - 3 - 8si(2si + 1)] 



(225) 



2(2s + l)(2s + 3) 



(226) 



c 2 o(g) + (4s 2 - 4s + 3)8si(2si + 1) 
2(4s 2 - l)(4s 2 -9) 



c 20 (s) = 3(4s 4 -9s 2 -s + 3), 



s[s 2 (4s 1 -3) + 8s 1 (s 1 + l-4s 2 )] 



4(4si - 1) 



sjJM + V) 
' (4^1-1) 



e?(s) = i [s 2 -4s 1 (2s 1 + l)] 



af(4ai + l) 
' (4si-l) 



(227) 
(228) 



(229) 



(230) 



for j — 1,2. For s% = 1/2, it is easy to see that 
aV 2 ' + (s) = by 2 ' + (s) = for s = 1, 2, as expected. 

Letting a/fis) = af, b' u± (8) = b ± , and c^f{s) = 
cf , it is easy to show that for Type I tetramers, 



1 



(231) 



3 

1 4si 

a„ + b~ = — a, = -Ci ,(232) 

2 3 1 3(4«i - 1) 1 ' V ' 

where = Cj x { . 

This implies that for Type I, the axial NN and NNN 
axial anisotropic exchange interactions may be combined 
to yield an effective axial anisotropic exchange interac- 
tion given by Eq. ([58]) . 

For the single-ion contributions to the level-crossing 
inductions, no such simple relation can be found. We 
note that 



a+ + & 4 



1 



(233) 



but the overall quantity a J + 2b + + cos 2 9 in Eq. (|202[) 
contains the extra quantity b + , which depends upon s, s±. 



H. Type II First-order AFM level-crossing constants 

For AFM tetramers with J g < 0, J' g - J g < 0, Type II, 
there are two classes of minimum energy configurations 
for each s manifold, depending upon whether s is even or 
odd. For even s, the relevant states have S13 = S24 = s/2, 
and for s odd, they are S13, S24 = (s±l)/2, (s=Fl)/2. This 
type is also relevant for FM tetramers with J g > J' g > 0, 
especially for the first excited manifold of states with 
s = 4si — 1. For even s the relevant parameters are 



Si, ± 
a s, s/2, s/2 

/i(s,si) 

"s, s/2, s/2 

/a («, si) 



1 



a -IT 



/i(s,si) 



2(2s- 1) L" * 1 2(s + 3) 
16si(si + 1) - s 2 - 2s + 6, 
1 



2(2s- 1) 



s 2 ± f 2(S, St) 



s + 3 

16si( Sl + l)(s 2 + 2s-l) 
-s(s 3 + 4s 2 + s - 4), 
s - 1 



-«,»/2,./a 2(2s-l)' 



(234) 
(235) 
(236) 

(237) 
(238) 



K>S! 

u s, s/2, s/2 



S 

64 



j(* + l)[16«i(*i + !)-«(* + 4)] 



16(s + 3) 
l)[16si(si + l) 



+ (1 "^ o) " 128s(s + 3) 

x(s + 2)[16s 1 (s 1 + 1) -s(s + 4)] 
4s 2 (s 1 + l) 2 

+Os,0 5 , 



4] 



s, s/2, s/2 



S 

64 



4 4s 2 ( Sl + l) 2 

Os,0 ^ 



(s+l)[16si(si + l) 



" (1 "^ o) 128s(s + 3) 

x(s + 2)[16si(si + l)-s(s + 4)]. 

For odd s, the relevant parameters are 



(239) 



4] 



(240) 



l s,(s+l)/2,(s-l)/2 



1 



2s(2s - 1) 



T: 



2(s + 2)(s + 4). 
/ 3 (s, Sl ) = 16s 1 (s 1 + l)(s 2 + 3s-l) 
-s 4 - 5s 3 + lis - 11, 
1 



\(j+l)/2,(»-l)/2 



2(2s- 1) 
± 



s 2 -l 



/i(s,si) 



(241) 
(242) 

(243) 



(s + 2)(s + 4)_ 
16si(si + l)(s 3 + 5s 2 + 4s - 3) 
-(s + l)(s 4 + 6s 3 + 7s 2 - 3s + 1), 



30 



+ _ s 2 - s + 1 

=«,(.+i)/a J (.-i)/a - 2s(2s- 1)' 

(s + l)(s-l) 2 
-«,(«+i)/a,{.-i)/a 2s 2 (2s-l) ' 

( s 2-l)2 ( s + i)2 



(244) 
(245) 



° s ,0+i)/a,(s-i)/a 



64 32(s + 2) 

x[16ai(»i + l)-(« + 3)(«-l)] 
+ [16s 1 (s 1 + l)-(s + 3)(s-l)] 2 

(s + l) 4 
* 256s 2 (s + 2) 2 

+[16ai(»i + l)-(a + 3)(«-l)] 
x[16«i(«i + 1) - (s + 5)(s + 1)] 

(2s + 3)(s 2 -2s-l) 2 
X 256s 3 (s + l)(s + 4) (2s + l) 2 ' 



(246) 



T-Sl 

s,(s+i)/a,(s-i)/a 



(, 2 -i) 2 ( S +ir 



64 256s 2 (s + 2) 2 

X[16s x (si + 1) - (s + 3)(s- l)] 2 
-[16si( Sl + l)-(s + 3)(s-l)] 
x[16si(si + 1) - (s + 5)(s + l)] 

(2s + 3)(s 2 -2s-l) 2 
X 256s 3 (s + l)(s + 4)(2s + l) 2 ' 

(247) 

From these expressions, we obtain the level-crossing in- 
ductions for the AFM type J g < and J' g - J g < 0. 
For s even, we have 

= s ( 2s - iK^/swa 

-( s -l)(2 s -3)Cj !/2 , M/2 , ( 248 ) 

si,± 

s,s/a,s/a 

(249) 



a IIe,2( S ) = Sa 



-(S l) a s ^i jS / 2 ,(s-2)/2' 



^00 = C/^-CUm/J' ( 25 °) 

c JIefl( S ) = *( 2a ~ ^^./^/a 

-(s-l)(2 S -3) C ±_ lis/2)(s _ 2)/2 , (251) 

C //o,2( S ) = SC ts/2,s/2 

-(s - l)c s _ l!S / 2) ( s _ 2 )/2' ( 252 ) 

^//e( S ) = ®7,s/2,s/2 ~ ^s-l,s/2,(s-2)/s' ( 253 ) 

e //e( S ) = ^7,l/2,s/2 _ •^-l,a/2,(*-2)/a- ( 254 ) 



Tie 

=>si 
'Jle 

For s odd, we have 



a //o,l( S ) — S ( 2s 1 ) a s!('s+l)/2,(s-l)/2 

-(s-l)(2s-3)a^ (s _ 1)/2(s _ 



l)/2> 



(255) 



Sl,± / \ Si.± 

a IIo,2\ S ) — Sa ;5,(s+l)/2,(s-l)/2 

-(«— 1 ) a ;5 L^( s _i)/2,( s _l)/2' ( 256 ) 



s,(s+l)/2,(s-l)/2 "s-l,(s-l)/2,(s-l)/2' 

(257) 

s(2s - l)C St ( s+ iy2,(s-l)/2 

— (s - l)(2s - 3)c s _ 1 ^ g _ 1 jy 2 ^ a _ 1 y 2 , 

(258) 



-7/o,2 



Jo, 



SC S ,(s+l)/2,(s-l)/2 

— (S - l)C s _ li ( s _ 1 )/ 2 ( s 

D s,(s+l)/2,( s -l)/2 

°s-l,(s-l)/2,(s-l)/2' 

T-Sl 

s,(s+l)/2,(s-l)/2 

'T'31 

s-l,(s-l)/2,(s-l)/2- 



-X)/2' 



(259) 
(260) 
(261) 



From these expressions, we may obtain the Type II 
AFM level-crossing induction parameters. For even s, 
we find 



a //M( s ) 
a c 10 (s) 



2 s~ 3 48s 1 (s 1 + l) + Q c 10 (s) 

2 T 4 (s + l)(s + 3) ' ^ j 
-2s 3 - 5s 2 + 6s + 18, 



(263) 



1 



2s 2 - 6s + 3 



2(2s-l)(2s-3) 
16si(si + l)(2s 2 - 4s + 3) + a| (s) 



± 



2(s + l)(s + 3) 



a 9n s 



J IIc 



(s) 



2s 4 + 2s 3 - 9s 2 - 18s + 18, 
s 



s - 1 



(2s- 1) (2s -3) 
16s x (ai + !)(«- 2) + &§(*) 



± 



TO 



'IZe,2 



2(s + l)(s + 3) 
-4s 4 - 7s 3 + 20s 2 + 14s - 18 
s(2s- 3) 

2 (*-l) ' 

s(2s 2 - 4s + 1) 
2(s- l)(2s- l)(2s-3)' 



(264) 
(265) 



(266) 
(267) 
(268) 

(269) 
(270) 



and d //c (s) and ej Je (s) are given by Eqs. (|253|) and (|254|) . 
Combining »//e~2( s ) an< ^ ^//c _ ( s )' we nn< ^ 



1 16si(si + l)-dg(s) 



, (271) 



"//c,2^//c - 2 4(s + l)(a + 3) 

d c {s) = 2s 3 + 7s 2 + 2s - 6. (272) 

We note that this expression differs substantially from 
that for a//e^i( s ) > except for si = 1/2 and s = 2. Simi- 
larly, it is elementary to combine cj lc 2 (s) and [b 6 //^ (s) + 
b%-(8)]/2. We find 



-7/o,2 



2(s-l) 



(273) 
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This simple expression differs from that for cj le x (s) by 
the factor 2s — 3. However, at s = 2, the only even s 
value for s\ — 1/2, they are equivalent. In addition, as 
for Type I, there is no simple relation between the single- 
ion parameters aj^'^is) + 2&jjj, + (s) and cijje~^(s). 

For odd s, the Type II level-crossing induction param- 
eters are 



I. Hartree INS functions 



l IIo 



\{s) 



a° (s) 



2s- 1 48si(s 1 + l) + a^ (s) 
2 T 4(s + 2)(s + 4) ' 
-2s 3 - lis 2 - 10s + 17, 
1 



2s z - 2s - 1 



2(2s- l)(2s-3) 
16si(si + l)(2s 2 + l) + a^ (s) 



± 



2(s + 2)(s + 4) 



a^ (s) = 2s 4 + 10s 3 + 9s 2 - 22s - 5, 

f (s-l)(s-2) 



(274) 
(275) 

(276) 
(277) 



J IIo 



is) 



(2s- l)(2s-3) 
16s 1 (s 1 + l)(s 2 -4s + l)+6g(s) 



± 



2(s + 2)(s + 4) 



TO 

"Jlo,l( S ) 
3 J7o,2( S ) 

dfio(l) 



-Ho 



(1) = 



-4s 5 - 19s 4 + 54s 2 - 2s - 5, 
(s-l)(2s-l) 
2s 

(s-l)(2s 2 -4s + 3) 
2s(2s- l)(2s-3) ' 

^(4s 3 + 8s 2 + 7s! + 3), 
6 

-^(4s 3 + 8s 2 + 3s!-l), 
6 



(278) 
(279) 

(280) 
(281) 
(282) 
(283) 



and the other djj (s) and ejj Q (s) values are found from 
Eqs. (|260j) and pTC]) . 



We note that a 



(,i/a, 



1/2,- 



(2) 



1/2,- 



(1) 



^' + (2) 



(1) = for j — 1,2, as expected. However, by 



combining etj^ 2 (s) and frj^ (s), we have 

1 16« x (ai + 1) - dg(a) 



a /7o,2 + fe //o" 



(284) 



2 4(s + 2)(s + 4) 
2s 3 + 13s 2 + 22s + 5, (285) 



which differs substantially from the version with si = 
1/2. In addition, it is elementary to combine cJ Io2 (s) + 

K/o + (s) + ^o>)]/2- We find 



'//o,2 



&-' + (s)+6— 



2s 



(286) 



which differs from cJ Io 1 (s) by the factor 2s — 1. At s = 1, 
the only relevant odd s value for s% = 1/2, these are 
equivalent. In addition, as for Type I and the even cross- 
ings of Type II, there is no simple relation between the 
single- ion parameters a^'^ls) + 2&^ + (s) and a//o~i(s). 



The functions L v y(q) and M u y(q) in the self- 
consistent Hartree INS Sg\q,u) in the induction rep- 
resentation are given by 



cr"=±l 



a"=±l 



S13 <-> S24 
s 13 ^ s 24 



(287) 



<T = ±1 

X 



er" = ±l 



2 Wh^'S 



',s+cr> \ U - (y ' s -(a' + l)/2 



o' = ±\ 



cr"=±l 



fv'M) 
fh(q) 

f±(q) 



Sl3 <-> S24 
s' 13 <-> S24 J • 

V % -» -% 

-2£ s , Sl3>S24 sin^a) sin^a) 
2(l - cos^fe + 

y ( pSl3 + (o-"+l)/2,S 2 4 
g/-('7) ? 7s + (cr' + l)/2,s 

2U - cos[a(9 ;c + q y )] 

( r is 13 + (a" + l)/2.s 2 4 
\y s 1 ,s 1 ,a'<j"s+a" (a' + l)/2) ' 

1 + cos ((fa a) cos(q y a) 
±cos(g z c)[cos(<faa) + cos(g y a)], 



(288) 

(289) 

(290) 
(291) 

(292) 
(293) 
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wlnerp tine A m C m V)° '> m P S 13,S24 r;s 13 ,s 2 4 „ £ 

are given by Eqs. g6]), (flBBf . (fl"67l) . and (fT72l) - (fT75| . 

respectively. 
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